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Abstract 



A systematic investigation of hard thermal photon spectra from central ultra-relativistic 
heavy ion collisions is presented with emphasis on the effects of bremsstrahlung processes 
in the quark-gluon plasma (QGP). Bremsstrahlung photon production in the quark-gluon 
plasma has recently been considered within the Braaten-Pisarski method in thermal QCD, 
where rates have been fomid that exhibit the same order in the coupling constants as those 
describing the lowest order processes, Compton scattering and gg-annihilation. The im- 
pact of these bremsstrahlung photon production rates on the thermal photon spectra is 
studied systematically within a simple, well understood one-fluid hydrodynamical model 
that describes an only longitudinally expanding fireball (1+1 Bjorken scaling hydrody- 
namics). A first-order phase transition is implemented in which QGP (simulated by an 
ideal massless parton gas of two-flavors) "hadronizes" according to the Gibbs criteria and 
Maxwell construction into a hot hadronic gas (HHG) (simulated by an ideal massless pion 
gas). It is found that the bremsstrahlung processes enhance the thermal photon yield 
from the QGP by about one order of magnitude over the complete considered p^-range 
independent of the choice of the model parameters. This results in an enhancement of 
the total thermal photon yield which is most significant for parameter sets that support 
a highly contributing QGP phase. The influence of each model parameter on the thermal 
photon spectra is examined carefully and a thorough understanding of the model is ob- 
tained. Experimental upper limits on direct photon production in fixed target 200 ^-GeV 
S + Au collisions at the CERN SPS are also considered and used to extract upper limits 
for the initial temperature of the QGP, where the QGP bremsstrahlung processes are 
found to make a difference of about 15 to 20 MeV depending on the temperature at which 
the phase transition is assumed. In comparison with other theoretical studies, the impor- 
tance of reaction features not described in the simple model are estimated and interesting 
elements for a future extension of this systematic investigation are identified, which will 
be of great interest in prospect of the upcoming experiments at the BNL RHIC and the 
CERN LHC. 
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Chapter 1 



Introduction 



This thesis presents a systematic investigation of thermal photons produced in ultra- 
relativistic heavy ion experiments. The primary goal of these experiments is the discovery 
of the quark-gluon plasma which is the dcconfined state of strongly interacting matter 
predicted by quantum chromodynamics. If nature allows the existence of this state, 
physicists face two challenges. First, they have to produce the quark-gluon plasma, and 
second, they must clearly identify that it indeed has been produced. For the first task 
high energy nucleus-nucleus collisions are the ideal means, for the second task many 
signals have been proposed. In this thesis the quality of thermal photons as a potential 
signature for the quark-gluon plasma will be examined theoretically. A systematic study 
of thermal photon emission will deliver insights into the space-time development of the 
fireball formed in ultra-relativistic heavy ion collisions. Within a simple hydro dynamical 
model for the fireball evolution the effects of the most recent thermal photon rates on 
the photon spectra will be analyzed. Of course, experimental photon data will also be 
inspected and compared with theoretical results. Due to decisive experiments underway 
the thermal photon investigation illustrated in this thesis will be of interest at least for 
the next decade. 

This chapter will provide the motivation for the systematic investigation of thermal 
photon production in ultra-relativistic nucleus-nucleus reactions, and it will give back- 
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ground information on heavy ion physics at ultra-relativistic energies. In Sec. 1.1, the 
accepted theory of strong interactions, quantum chromodynamics, will be reviewed. This 
theory predicts quark-gluon plasma as the state of nuclear matter at high temperatures 



and high densities. Section 1.2 will describe this state of matter and its presumed appear- 
ance in nature. In the subsequent section, Sec. 1.3, the pursued production of quark-gluon 
plasma in the laboratory will be discussed. Then Sec. |l.4| will center on proposed signa- 
tures for the quark-gluon plasma. Because the focus of this thesis is on photons, the 
emphasis will be on distinguishing properties of electromagnetical probes. The experi- 



mental situation, present and future perspectives, will be the topic of Sec. |1.5| . Based on 
the background information gathered in these preceding sections, the motivation for the 



systematic study of thermal photons will be addressed in Sec. 1.6. An overview of this 
thesis will conclude Chap. ||. 



1.1 Quantum Chromodynamics 



In our present understanding of nature, quantum chromodynamics [QCD) is the theory 
of strong interactions [j^. It is a non-Abelian gauge theory which is based on the SU(3) 
color gauge group. Thus, the fundamental principle is local gauge invariance under SU(3)- 
transformations. For keeping up this symmetry, gauge fields are crucial. The gauge quanta 
of QCD are the gluons, massless bosons of spin one. Gluons are exchanged between 
particles that carry color, the quantum number of QCD. Such gluon sources are the 
quarks, massive particles of spin one-half. Today, six different quark types or "flavors" 
have been observed, the up (u), charm (c), and top (t) quarks having electrical charge 
+2e/3, and the down (d), strange (s), and beauty (b) quarks having electrical charge 
— e/3. Every quark, independent of its flavor, comes in one of three colors, e.g., red, blue, 
or green. However, not only quarks can emit and absorb gluons, the gluons themselves 
carry one color and one anti-color. Because a gluon cannot be color-neutral, there are 
as many as eight gluons with different color "charge." Consequently, the gauge fields of 
QCD interact among themselves. It is this distinguishing property that made QCD the 
accepted theory of strong interactions. 
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There are two phenomena characteristic for the strong force: asymptotic freedom and 
confinement. High energy deep inelastic scattering experiments revealed asymptotic free- 
dom: at close distances, quarks behave like free particles. QCD can derive this observation 
from first principles. As long as there are no more than 16 flavors, QCD is asymptotically 
free: the coupling constant becomes weak at high energy. This decrease of the strong 
coupling constant at high energy implies an increase at low energy. It is a hypothesis 
inferred from these considerations, that quarks or gluons cannot be observed as isolated 
particles. This hypothesis, called confinement, matches the experimental fact that nei- 
ther an isolated quark nor an isolated gluon have ever been detected. Only color-neutral 
particles can be found as isolated objects in our physical world. 

Baryons and mesons are color-neutral particles made up of quarks and gluons. In a 
baryon, the combination of three constituent quarks, each with a different color, leads to 
a color-neutral state. In a meson, the color-neutral state is realized in the combination 
of a colored quark and an anti-quark bearing the corresponding anti-color. Of course, 
one can think of various other combinations that lead to color-neutral objects, such as 
glueballs or mesonic molecules. These exotics have not been detected clearly, however, 
their observation would confirm QCD as the theory of strong interactions, a fact that 
motivates an area of current research in nuclear physics ||2|. 

1.2 The Quark-Gluon Plasma 

There is another exotic state of matter predicted by QCD, the quark-gluon plasma (QGP) 
iP, ^, ^, p. According to QCD, the strong coupling constant decreases at high energy. 
Thus, in an environment of extremely high energy density, quarks and gluons are expected 
to form a relativistic weakly-interacting parton gas: a QGP. The search for QGP is one 
of the central topics in strong interaction physics, and also this thesis aims at reviewing 
the QGP formation in heavy ion collisions. 

The exploration of the QGP gains much attention because it is considered a key to 
various fundamental questions. For example, the standard model of cosmology assumes 
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the existence of a QGP phase in the early universe. It is beheved that this color deconfined 
state of the universe underwent a phase transition about 10~^ seconds after the big-bang, 
where the quarks and gluons became confined in baryons and mesons. This process, called 
hadronization, is directly connected to confinement and far from being understood. It is 
manifested in high-energy physics experiments as jet formation, but can only be repro- 
duced in simplified theoretical models 0. In the universe we see today, QGP could still be 
present, since supernovae and neutron stars provide extreme astrophysical environments 
which favor the creation and existence of QGP |8|. 

Presently, nuclear matter under extreme conditions can be produced and studied in 
heavy ion collisions [l^. These experiments are ideal means to compress and heat up 
nuclear matter in the laboratory, where they give a unique tool to determine the nuclear 
equation of state. If the QCD predictions are valid, this nuclear equation of state should 
contain a deconfinement phase transition to QGP at high temperatures and densities as is 
illustrated in Fig. In fact, physicists are confident to produce QGP in ultra-relativistic 
heavy ion collisions if nature allows its existence. 



1.3 Ultra- Relativistic Heavy Ion Collisions 



Ultra-relativistic heavy ion collisions are performed at the Brookhaven National Labora- 
tory (BNL) and at the European Center for Nuclear Research (CERN), where heavy ions 
are accelerated in high-energy proton accelerators before they are aimed at fixed targets. 
In this way the BNL Alternating Gradient Synchrotron (AGS) has provided "^^Si —^^'^Au 
and ^^"^Au —^^"^Au collisions with a center-of-mass energy of up to ^/s = 5 A ■ GeV. At 
CERN, the Super Proton Synchrotron (SPS) is still producing ultra-relativistic heavy ion 
collisions in fixed-target experiments, center-of-mass energies of up to ^/s = 20 A ■ GeV 
are reached for —'^^^Au and ^^^Pb —^^^Pb reactions. In the experiments at AGS and 
the medium-energy experiments at SPS, a complete stopping of the baryonic projectile 
constituents in the middle of the reaction zone was observed. The production of a short- 
lived QGP is expected because the significant overlap of the baryons presumably causes a 
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pc p 



Figure 1.1: Phase Diagram of Nuclear Matter. At liigh temperatures T and high densi- 
ties p, nuclear matter is expected to be in the QGP state. By lowering the temperature 
and the density, one shoTild cross the critical boundary where quarks and gluons become 
confined into hadrons. Below the phase transition values of temperature, Tc, and density, 
Pc-, we experience nuclear matter in its hadronic state. 
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screening of the color confining potential. Experiments in this energy region, the stopping 
region, are considered ideal for delivering insights into the nature of neutron and hybrid 
stars for which astrophysicists predict similar baryon dense conditions. A schematical 
drawing of a central heavy ion collision in the stopping region is shown in Fig. |1.2| (a). 



A different energy region, the transparent region, is possibly already considered in the 
high-energy nucleus-nucleus collisions at SPS, and it will definitely be studied in the next 
generation of ultra-relativistic heavy ion experiments. By applying the collider principle, 
the Relativistic Heavy Ion Colhder (RHIC) at BNL will achieve ^ = 200 ^ • GeV with 
^^"^Au projectiles. At the moment, this dedicated heavy ion accelerator is in its testing 
phase. In the year 2005, even higher energies will be available with the completion of 
the Large Hadron Collider (LHC) at CERN. This device will collide ^^^Pb projectiles up 
to center-of-mass energies of ^/s = 5500 A • GeV. Because of the high energies pursued 
in these upcoming experiments, the accelerated heavy ions will suffer extreme Lorentz 
contraction, therefore, the overlap of the projectiles will not last long enough to stop the 
nuclei significantly. The heavy ions will instead be approximately transparent keeping 
much of their initial energy 12]. However, there will be a strong color field between 



the emerging baryons of the initial projectiles, which will polarize the vacuum and cause 
parton pair production. In this transparent energy region, physicists expect the produc- 
tion of a quasi baryon free QGP. As the QGP phase of our early universe was also baryon 
free, heavy ion reactions at these energies should help in understanding the QGP phase 
of our early universe. Figure |1.2| (b) illustrates a central ultra-relativistic nucleus-nucleus 
collision in the transparent region. 

This thesis concentrates only on transparent heavy ion reactions because the central 
region of the longitudinally expanding fireball is ideally suited for the study of thermal 
photon yields. A significant stopping of the baryons of the initial projectiles and conse- 
quently a finite baryon density in the central region is already expected in SPS high-energy 
nucleus-nucleus collisions. However, for simplicity, we also treat high-energy SPS heavy 
ion experiments as being in the transparent region and keep in mind that this might only 
be a fair approximation. 
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(a) 





(b) 



I 






! 



Figure 1.2: Ultra-Relativistic Heavy Ion Collisions. According to the center-of-mass en- 
ergy ^/s, one distinguishes the stopping (a) and the transparent region (b). The two 
scenarios shown are, of course, only the limiting cases. For example, in high-energy SPS 
heavy ion experiments a scenario in between these pictured seems realistic. 
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A common picture for the space-time evolution of a transparent ultra-relativistic heavy 



ion colUsion j^S] is illustrated in Fig. 1.2. The maximum overlap of the colliding nuclei, 
which are highly Lorentz contracted, defines proper time r = fm in the considered 
system. By multiple scatterings among the initial partons produced through vacuum 
polarization, the system goes into thermal equilibrium at initial time tq. In the transparent 
energy region, the initial temperature Tq is expected to be sufficiently high for supporting 
the existence of a QGP. As the nuclei emerge, the system expands and cools. At time 
the transition temperature Tc is reached, where the phase transition sets in. Assuming 
a first-order phase transition, a mixed phase (MP) of constant temperature Tc follows, 
in which quarks and gluons become confined into hadrons. The system cools further 
when all the quark matter has transformed into hadronic matter. This completion of 
the phase transition takes place at r^. The produced hadrons do still interact among 
themselves forming a hot hadronic gas (HHG). When the system arrives at the freeze-out 
temperature Tf, the hadrons will stream as free particles out of the collision zone. 



1.4 Signatures of the Quark-Gluon Plasma 

Provided nature allows the existence of QGP, the challenge of ultra-relativistic heavy 
ion experiments will not only be the production but also the clean identification of this 
deconfined state. An ideal signature would be an observable physical phenomenon that 
can only be explained by assuming the production of a QGP phase. This signature 
should not be understandable in any realistic model not embedding the deconfinement 
phase transition. Thus, potential signals need to be examined in both phase transition 
and no phase transition scenarios. 

Many phenomena have been proposed as signals for QGP formation [p!^]. They can 
be grouped into hadronic and electromagnetic observables. Hadronic signatures are, e.g., 
strangeness enhancement, J/ip -suppression, and detection of strangelets, which are ex- 
otic objects that contain several strange quarks. Also, thermodynamic signals can be 
considered hadronic because thermodynamic variables, such as temperature T, entropy 



8 



Figure 1.3: Space-Time Evolution of an Ultra-Relativistic Heavy Ion Collision in the 
Transparent Region. The Minkowski-diagram in the right half of the figure drafts the 
fireball behavior along the beam-axis labeled by z as seen in center-of-mass time t. The 
left half of the figure sketches the spatial evolution of the heavy ion reaction on the 
same time-scale. The projectiles are almost light-like moving practically on the light- 
cone. At the point of maximum overlap, (t, 2) = (0,0), the initial nuclei are hardly 
stopped. Then, a strongly interacting continuum is produced by the strong color field 
of the receding projectile fragments. The hyperbolas in the upper time-like part of the 
diagram characterize domains of constant proper time. These domains mark the starting 
and ending points of the different collision phases as stated in the text. 
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density s, and energy density e, are inferred respectively from average transverse momen- 
tum, multiplicity and transverse energy distributions of hadrons. In general, hadronic 
observables can be measured easily due to the copious production of hadrons in ultra- 
relativistic heavy ion reactions. However, with the exception of jets, hadrons cannot 
carry direct information about the early collision stage because they interact until the 
system undergoes freeze-out, as mentioned in the previous section. Therefore, observed 
hadrons carry primarily information on the freeze-out stage of the collision. For elec- 



tromagnetic observables, photons and dileptons, the situation looks different |15|. Since 
electromagnetic mean free paths are much larger than the transverse size of the fireball, 
most photons and dileptons produced in the reaction reach the detector with no final state 
interaction. In this way even the earliest and hottest phase of the collision can be probed. 
However, in experiments electromagnetic signals are hard to detect. The overall rate of 
electromagnetic probes is small and needs to be extracted out of large backgrounds from 
hadronic decay processes. Nevertheless, we concentrate on electromagnetic signals: em- 
ploying a well-understood model for the space-time evolution of the fireball, we investigate 
systematically the production of thermal photons that have momenta above the fireball 
temperature, > T. 



1.5 The Experimental Situation 

Most important within a systematic theoretical study is a comparison with experimental 
results. For thermal photon investigations experiments provide direct photon yields. Di- 
rect photons are those produced inside the strongly interacting continuum present before 
freeze-out. In their extraction, a clean separation of decay photons is crucial. These 
photons mainly originating from Dalitz decays vr^ — > 77 and 7/ ^ 77 dominate the ob- 
served photon spectra. The measurement of direct photon spectra was a main goal in 
the experiment WA80 at the CERN SPS. In fact, the only reliable direct photon data 



has been published by the WA80 collaboration |16, 17 1, where upper limits on the direct 
photon production in 200 A ■ GeV S + Au collisions were determined on a statistical ba- 
sis. At present, the WA80 experiment has been upgraded to WA98, which is examining 
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158 A ■ GeV Pb + Pb collisions at the CERN SPS. Because of the larger projectiles higher 
yields not only of direct photons but unfortunately also of resonance decay photons are 
expected. However, with the sophisticated WA98 photon spectrometer LEDA (LEadglass 
Detector Array) the WA98 results for direct photon production should be more precise 
than those from WA80. Today, the results of the WA98 direct photon analysis are eagerly 
awaited and will be published in the very near future. For results from RHIC and LHC, 
we have to be patient. PHENIX, one of the two large-scale detectors at RHIC, will start 
data taking in fall 1999. With high granularity and excellent particle identification capa- 
bilities, it will be well suited to extract the direct photon signal in the high-multiplicity 
environment predicted for the RHIC experiments [^], but it might take about five years 
until the PHENIX collaboration can provide more than preliminary direct photon pro- 
duction data. Then, in 2005, the CERN LHC will start operating. At this accelerator, 
there will be ALICE (A Large Ion Collider Experiment), a dedicated high energy heavy 
ion physics experiment. ALICE will be equipped with PHOS, a state-of-the-art electro- 
magnetic calorimeter designed for photon physics in high multiplicity reactions ||l^ . With 
the experience collected in WA80, WA98, and PHENIX, ALICE will give presumably the 
cleanest extraction of direct photons, yet we will need to wait another decade for these 
results. 

1.6 Motivation for a Systematic Investigation of Thermal 
Photons 

Our research ultimately aims at the discovery of the QGP and the connected confirma- 
tion of QCD by considering hard thermal photon production in ultra-relativistic heavy 
ion collisions. For a theoretical prediction of the photon spectra measured experimen- 
tally, cross sections of the elementary photon emitting reactions, or photon production 
rates, need to be calculated and integrated over the space-time history of the fireball. 
Recently, thermal photon rates have been revised for the QGP phase in the framework of 
finite temperature QCD [|2^. An astonishing result has been found: because of medium 
effects, bremsstrahlung processes contribute in the same order of the coupling constant as 
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the Compton scattering and gg-annihilation processes calculated in ||2l|, This leads 
to a significantly enhanced photon production rate. In fact, bremsstrahlung processes 
become the dominant source for hard thermal photons in the QGP phase. Earlier sys- 
tematic investigations implemented for the QGP state only the Compton scattering and 
gg-annihilation rates |23, 24, 25, which means thermal photon emission of the QGP 
phase was underestimated. A further investigation with emphasis on the effects of the 
new rates seemed very necessary, and it was this consideration that triggered the research 
for this thesis. Because experiments cannot distinguish in which stage of the fireball a 
detected photon was emitted, thermal photons from the HHG phase were also included in 
our study. We extract these photon spectra employing the most recent parameterization 
of the production rate in hadronic matter p7| , |2^ . For the dynamics of the fireball, many 
models can be found and more are under construction. The variety reflects the ongoing 
discussions and uncertainties on the space-time evolution of an ultra-relativistic heavy 
ion collision. Since our prime interest is on the effect of the rates due to bremsstrahlung 
processes in comparison to the rates used before (Compton scattering, gg-annihilation) , 
we employ only a simple hydrodynamic model. The physics basis of this model is well 
understood and the effects of the rates will not be covered by fancy features of the reac- 
tion dynamics. However, the price for good understanding are the severe limitations one 
has to bear in mind, e.g., the model will not describe any transverse expansion or any 
non-equilibrium behavior. Within this model that has a decent number of parameters, 
we perform the systematic study on thermal photon spectra. Direct constraints on the 
parameter set will be inferred from WA80 photon data. More severe constraints are ex- 
pected from comparisons with future experimental data of WA98, PHENIX, and ALICE. 
With these upcoming constraints we are confident that our systematic investigation will 
contribute in the verification of the QGP state. 



1.7 Overview of Thesis 



The next chapter will describe the simple model we use for the space-time evolution of the 
fireball. A short review on relativistic hydrodynamics and Bjorken initial conditions will 
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be given followed by a discussion of the nuclear matter equation of state. Chapter |3| will 
be dedicated to photon production in hot thermalized, strongly interacting environments. 
We will enumerate and investigate photon production processes embedded in our calcu- 
lations. Having studied the current thermal photon rates, the results of our systematic 
investigation will be presented in Chap. ^. This part will address separately the influence 
of each model parameter. Comparisons with other works and results from the analysis of 
experimental data on direct photon production will be given in Chap. |5[ 
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Chapter 2 

A Simple Model for 
Ultra-Relativistic Heavy Ion 
Collisions 



For a systematic investigation of thermal photons as a potential signature of QGP, one 
must convolute the rates with the space-time history of the nucleus- nucleus collision. The 
elementary photon rate, E^^ dN / {d!^x d?p^,), is the number of photons dN with energy 
emitted per unit volume per unit time within the three-momentum interval [p-y, p^+d^p^]. 
For thermal photons, this rate depends, of course, on the temperature T of the emitting 
space-time point, 

^ dN ^ dN 

Since photons cannot be traced back experimentally to their origin in space-time, only 
energy and three-momentum are measured. Consequently, for a comparison with the 
detected photon spectrum, the rates need to be integrated over the space-time evolution 
of the heavy ion reaction. 

The integration limits and the temperature field T{x) are necessary for the evaluation 
of this four-dimensional integral. While some integration limits are determined by the 
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accelerator and the projectile nuclei, others can be considered parameters. However, the 
temperature field T{x) should be delivered as a feature of some reaction model. We obtain 
this field within a simple hydro dynamical picture of an ultra-relativistic heavy ion collision 
that is described in this chapter. 

We use relativistic one-fluid hydrodynamics because it provides a simple but, in our 
view, realistic framework to study thermal photon production in high-energy nucleus- 
nucleus reactions. It is our strategy to employ a model that can be well understood and 
maybe even solved analytically. Within this approach, we find very clearly the influences 
of the photon rates and the model parameters on the thermal photon spectrum, which 
will be presented in Chap. 4. The influence of the nuclear equation of state (EOS) can 
also be checked because it is a crucial input of a hydrodynamical calculation. 

Our idea of a relativistic heavy ion reaction embedding a deconfinement phase tran- 
sition was already discussed in Sec. |1.3[ It is this picture we refer to whenever we talk 
of a phase transition scenario. On the presented space-time development, the use of a 
hydrodynamical simulation seems reasonable. A stringent condition for a system to be 
described by hydrodynamics is local thermal equilibrium. This allows a hydrodynamic 
simulation of the fireball evolution only between initial time tq and freeze-out time rj. 
Non-equilibrium behavior present before tq, also called thermalization time, cannot be 
modeled. In addition, the phase after freeze-out in which the particles are heading freely 
towards the detector is not suited for a hydrodynamical description. However, because 
the emphasis is on thermal photon spectra, no thermal photon yields are lost from these 
restrictions. Thermal photons can, of course, only be emitted from a thermalized system. 

In the subsequent section, we introduce relativistic hydrodynamics. To perform a 
hydrodynamical calculation, important ingredients are necessary as the EOS and initial 
conditions. In the simulation, we use Bjorken initial conditions, which allow a scaling 
ansatz of the four- velocity field u^{x) that describes the hydrodynamic flow. Together 
with a simple EOS, this scaling ansatz enables an analytical derivation of the temperature 
field T{x). While Bjorken initial conditions are covered in Sec. |2.2| , the EOS will be 



addressed in Sec. 2.3. We complete the description of the employed hydrodynamical 
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model by discussing two scenarios, one with and the other without a deconfinement phase 
transition. 



2.1 Relativistic Hydrodynamics 



This section discusses the essentials of relativistic hydrodynamics [29|. Within relativistic 
hydrodynamics, strongly interacting matter is considered a relativistic fluid. Because we 
neglect any dissipative effects as heat transfer or viscosity, it is even treated as a perfect 
or ideal relativistic fluid. The use of relativity is required because the particles in the 
collision system have velocities close to the velocity of light. There are two sorts of ve- 
locities considered in hydrodynamical calculations; collective (macroscopic) velocities and 
thermal (microscopic) velocities. Thermodynamic functions, such as energy density e and 
hydrostatic pressure P, characterize the thermal motion of the microscopic constituents 
of the continuum. Thus, they are a measure for the thermal velocities. The collective 
velocities describe the hydrodynamic flow in the form of the four-velocity field 

u^' = ^il,v) (2.3) 

with 

(2.4) 



V 



2 



where v is the spatial flow velocity vector. By checking 

u^W = 1, (2.5) 
one can easily see that is a time-like unit vector. We follow Landau's approach in 



attaching the hydrodynamic flow to the flow of energy [30|. In other words, is always 
tangential on the world lines of the energy flow. This defines the local rest frame (LR) 
of the fluid as the frame in which the energy flux vanishes. The collective velocity in this 
frame has the form 

nf^^) = (1,0,0,0). (2.6) 
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The hydrodynamic equations of motion are basic conservation laws. One is local 
conservation of energy-momentum 

9,,r^^(x) = 0. (2.7) 

In the above expression, T^^'^ denotes the energy-momentum tensor, which can be written 
for perfect fluids as 

r^^(x) = [e{x) + P{x)] u^{x)u''{x) - P{x)g^'^ (2.8) 

with the metric tensor g'^'^ = diag(l, —1, —1, —1). Further, for any conserved scalar quan- 
tity, i.e. baryon number, a continuity equation holds locally for the corresponding volume 
density p, 

df, [p{x)u''{x)] = 0. (2.9) 



Another important conservation law can be derived by contracting Eq. (2.7) with Ui,{x). 
Applying the first law of thermodynamics on the contracted equation exhibits local con- 
servation of entropy 

[s{x)u^'{x)\ = 0, (2.10) 

where s denotes entropy density. 



In order to solve the equations of motion, ( |2.7] ) and (|2.9|), an additional equation of the 
form P = P(e,p), the EOS, is needed. Together with the EOS, the equations of motion 
form a closed system. However, before solving, initial conditions must be specified. 



2.2 Bjorken Initial Conditions - The Bjorken Model 

In this thesis, energy regions are investigated in which nuclear transparency is present. 
One-fluid hydrodynamics cannot describe this "leading baryon" effect, and two-fiuid |31] 



and three-fiuid [32| models are too technical to fit our strategy. Therefore, sticking to 
the simple one-fiuid model, the transparent behavior of the colliding nuclei has to be 
implemented in the initial conditions. Concentrating on the central region in which the 
baryon density is small and can be neglected and ignoring the fragmentation regions that 
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contain most of the baryons close to initial velocity greatly facilitates the solution of the 
hydrodynamic equations. 



Following Bjorken's approach |11], we describe only the longitudinal expansion of sys- 
tems produced in central ultra-relativistic heavy ion collisions. Any transverse expansion 
is neglected! Thus, we employ a 1 + 1 dimensional model considering an only longitudi- 
nally expanding tube of strongly interacting matter. Within this model the tube radius 
appears as an integration limit in Eq. (p.2|). Because only central collisions with zero 
impact parameter are studied, we take for this radius the projectile radius as given by the 
simple phenomenological formula 

Ra = 1.3 fm (2.11) 



Because "fluid" cells within the longitudinally expanding tube move with relativistic 
velocities, it is sensible to describe the fluid cells in the variables, proper time r, and 
rapidity y. For an expansion only parallel to the beam-direction, these variables have the 
following form, 

r = Vt^ - (2.12) 

1 , t + z , , 

y = -,'o,—, (2.13) 

where t and z are the center-of-mass frame coordinates for time and longitudinal po- 
sition of the corresponding fluid cell, respectively. In this reference frame, the origin, 
(t, X, y, z) = (0, 0, 0, 0), was chosen as the point where the two incident nuclei have maxi- 
mum overlap. We showed already in Fig. curves of constant proper time as hyperbolas 
in a Minkowski-diagram illustrating the space-time evolution of a heavy ion collision. The 
rapidity y is the variable that indicates the position of the fluid cells on these hyperbolas. 
The proper time r is just the variable coinciding with the local time t(L/?) iii the rest 
frame of the considered fluid element. 

Of course, the initial conditions should also be given in these variables. Motivated by 
the observation of a "central-plateau" structure in nucleon-nucleon and nucleon-nucleus 
collisions, Bjorken asserted the initial conditions being invariant under Lorentz boosts 
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along the beam-axis. We follow this assertion which means that the initial conditions 
imposed at initial proper time tq do not depend on rapidity y. Since the hydrodynamic 
equations are Lorentz covariant, this symmetry is preserved for the complete hydrodynam- 
ical evolution. More specific, thermodynamical quantities, as energy density e, pressure P, 
temperature T or entropy density s, remain independent of rapidity y, 

e = e(T), P = P{t), T = T{t), s = s{t). (2.14) 

In this picture of an only longitudinal expansion exhibiting the symmetry discussed 
above, a scaling ansatz for the four-velocity can be made, 

u''{t,z) = -{t, 0,0, z) . (2.15) 
r 



With this ansatz and conservation of energy- momentum (^J) , the basic differential equa- 
tion of Bjorken's hydro dynamical model, 

de e + P , , 

— + = 0, 2.16 

or r 



can be derived. Further, by inserting ( 2.15| ) into ( 2.1C ) the equation for entropy conser- 
vation gets the form 

? + ^ = 0, (2.17) 
dr T 

or equivalently 

^{st) = 0. (2.18) 

The two latter equations imply that the entropy per rapidity slice dS/dy is constant of 
the motion [11|, 

which can be used to infer the initial entropy density, sq = s{tq), from the measured 
multiplicity distribution dN/dy. 



2.3 The Equation of State 

As emphasized in the preceding sections, the EOS is an important ingredient in any 
hydrodynamic calculation. Together with the EOS the hydrodynamic equations become 
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deterministic. For strongly interacting matter, the EOS should be derived directly from 
the QCD Lagrangian. Because of the nonperturbative character of QCD at large spatial 
distances, this task is far from being trivial and an analytical derivation of the EOS from 
first principles seems impossible. Instead, one performs computer simulations of QCD 
on a discrete lattice of space and time. These lattice gauge calculations, which are built 
on first principles, work in the nonperturbative regime of QCD. The EOS of strongly 



interacting matter can thus be inferred quantitatively from lattice QCD 34 1. In fact, 
this numerical approach provided a first quantitative prediction of the deconfinement 
phase transition [^]. However, we employ only a model EOS that displays features of 
lattice QCD results. We construct EOS's separately for both deconfined and confined 
matter. Assuming a first-order phase transition, we match the two EOS's by Maxwell 
construction. Within this approach of an idealized EOS, we again follow the philosophy 
of keeping the physics basis of our simulation well understood. 

Because only the central region of ultra-relativistic heavy ion collisions is investigated, 
we neglect baryon density Pb{x) completely and assume chemical equilibrium. By con- 
sidering an only longitudinally expanding fireball, we are left with one hydro dynamical 
equation which is the basic differential equation of Bjorken's hydrodynamical model ( p. 16 ). 



Consequently, for solving the hydrodynamical calculation, an EOS of the form P = -P(e), 
or equivalently e = e{P), needs to be specified. A bag model EOS for the QGP phase and 
an equally simple EOS for the hot hadronic gas (HHG) phase is presented. Finally, the 
issue of the phase transition, which we implement as a first-order transition, is addressed. 



2.3.1 Quark-GIuon Plasma - The Ideal Massless Parton Gas 

A prerequisite for the formation of a QGP is an extremely high energy density, e > 
1 GeV/fm^. In this environment, the QCD coupling constant tends to zero, a phe- 
nomenon already mentioned as asymptotic freedom, and quarks and gluons form to a 
good approximation a noninteracting relativistic quantum gas. This can be seen in lattice 
QCD calculations, where the energy density follows the Stefan-Boltzmann law for tem- 
peratures higher than about 2Tc. Because in a QGP, light quarks, as u, d, and maybe 
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s, would dominate, quark masses can be neglected and the description of the QGP as an 
ideal massless parton gas seems reasonable. In calculating the thermodynamic quantities 
of this ideal massless parton gas, one needs to regard the different statistics that quarks 
and gluons obey. While Fermi-Dirac statistics governs the thermodynamics of quarks, 
gluons follow Bose-Einstein statistics. The nature of the vacuum in which an ideal parton 
gas can exist must also be taken into account. This is simplest done by giving this pertur- 
bative QCD vacuum a constant energy density B, known as bag constant. Typical values 
lie in the region around B'^/'^ = 200 MeV. The bag constant appears also with a different 
sign in the pressure at the boundary of the QGP, where it describes phenomenologically 
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confinement of the partons within the QGP "bubble" |36]. Under the above considera- 
tions, the following Stefan-Boltzmann expressions can be derived for zero quark chemical 
potential, ^q = 0, 

(2.20) 
(2.21) 
(2.22) 

where QGP as the referred state of matter is indicated in the subscript q. gq is the effective 
number of degrees of freedom. For QGP with Nc colors and Nf flavors, 

gq = 2[Nl-l) + {^-^AN,Nf (2.23) 

and with the standard value Nc = 3, one gets gq = 37 for a two-flavored and gq = 47.5 
for a three- flavored QGP. Finally, from the expressions for pressure ( 2.2(]| ) and energy 
density ( 2.21| ), one can directly read off the bag model EOS 



Eq = 3Pq + iB. (2.24) 



2.3.2 Hot Hadronic Matter - The Ideal Massless Pion Gas 

While in the limit of high temperatures strongly interacting matter can be considered an 
ideal massless parton gas, an equally simple description can be realized for low temper- 
atures. The picture of hadronic matter as an ideal massless pion gas should be a decent 
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approximation at temperatures Tc > T > rriT^. Taking into account finite hadron masses 
and also more massive hadrons than pions would definitively give a HHG EOS that is 
closer to reality. However, since finite mass states are Boltzmann-suppressed, the higher 
number of hadrons considered would be compensated to some extend. Due to this fact and 
the numerical demand necessary to get the EOS in this more realistic picture, we apply 
the ideal massless pion gas where the EOS can be derived analytically. This system gov- 
erned by Bose-Einstein statistics, of course, also exhibits Stefan-Boltzmann expressions 
for the thermodynamic quantities, 

Ph = (2.25) 

Sh = 9h^T\ (2.26) 
2vr2 

Sh = 9h^T\ (2.27) 

where the index h specifies the considered state of matter as HHG and Qh denotes the 
effective number of degrees of freedom. For a hadronic gas consisting of only pions, 

9h = 3. (2.28) 

The equation of state embedded in Eq. ( ^.25 ) and ( |2.26| ) is the well known ideal gas EOS 
for massless particles, 

eh = 3Ph. (2.29) 

2.3.3 The Phase Transition 

The order of the deconfinement phase transition is still a matter of ongoing research 1 33 ] . 



However, we implement a first-order phase transition in our "phase transition scenario" 
that implies a mixed phase in which QGP and HHG coexist. During this mixed phase 
the ideal massless parton gas "hadronizes" continuously into the ideal massless pion gas 
at a hadronization rate that can be characterized in the volume fraction of the QGP, 

Mr) - - (2 30) 

where Vq and Vh are the spatial volumes occupied by QGP and HHG, respectively, and 
Vtot denotes the total spatial volume of the fireball. The quantity A equals one during 
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the pure QGP phase and at the onset of the hadronization, tq < t < r^, and equals zero 
during the pure HHG phase and at the offset of the hadronization, < r < r/ . 

The thermodynamic properties at the critical boundary are determined by the Gibbs 
criteria, 

T§ = T^} = Tc, (2.31) 
P,^ = = Pc, (2.32) 

which express thermal and mechanical equilibrium, respectively. The superscripts mark 
the state of matter whose thermodynamic quantity is considered; q indicates the QGP 
state and h the HHG state. With the above equations and the expressions (^ ) and (|2]2|), 



the critical or transition temperature can be obtained in terms of the bag constant B and 
the effective degrees of freedom, and gh, 



Tc={ L (2.33) 



Figure 2.1 illustrates this relation for gq = 37 (two-flavored QGP) and gh = 3 (ideal 
massless pion gas). For example, using the reasonable value i?^/^ = 200 MeV, one obtains 
Tc = 144 MeV. In a systematic study one should vary the bag constant i3 as a funda- 
mental parameter, however, we will vary the parameter Tc, which is equivalent because 
of Eq. ( ^.33 ). While temperature and pressure remain constant at the critical boundary. 



energy and entropy density decrease. The evolution of the entropy density is given simply 
by Eq. ( 2.1^ ), which can be restated as 



sir) = (2.34) 
r 

In fact, this equation determines the time scale of the complete hydrodynamic expansion. 
A different expression for the entropy density in the mixed phase Sc can be obtained by 
Maxwell construction, 

s,{t) = A(t) sqir!) + [1 - A(r)] Shir^), (2.35) 



which can be used together with Eqs. ( |2.22| ), ( |2.27 ), and (2.34) in deriving an explicit 
form of the QGP volume fraction 
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Figure 2.1: The Dependence of the Critical Temperature on the Bag Constant. Expres- 
sion ( p.33| ) is shown graphically for Qq = 37 (two-flavored QGP) and gh = (ideal massless 
pion gas). In Chap. ^ we will treat Tc as a parameter. Because Tc is in fact determined 
by B, we only imply the more fundamental setting of B. 

It can easily be checked that this equation reproduces the properties of A discussed above. 
The proper time dependence of A is illustrated for different numbers of effective degrees 
of freedom in Fig. |2.2| . As a consequence of entropy conservation, the hadronization rate 
rises with a rising number of effective hadronic degrees of freedom g^. Conservation of 
entropy also governs through A the evolution of the critical energy density Ec, for which 
the expression 

e,(r) = A(r)e,(r,«) + [1 - A(r)] suir^) (2.37) 
is obtained by performing again a Maxwell construction. 

2.4 The Phase Transition Scenario 

In the preceding sections, we illustrated qualitatively our idea of an ultra-relativistic heavy 
ion collision. Since we extract thermal photon yields from a model which is based on this 
picture, we now present a quantitative discussion. 

25 
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Figure 2.2: Evolution of the QGP Volume Fraction. The proper time dependence of the 
QGP volume fraction A is illustrated for = 3, 9, and 27. The only other parameter 
important for the shape of A is the effective number of degrees of freedom present in 
the QGP, which was set to the value of two-flavored quark matter, Qq = 37. Because 
we consider an adiabatic expansion, the entropy contained in the QGP must be fully 
transferred to the constituents of the HHG. For a high number of effective hadronic 
degrees of freedom, e.g., gh = 27, which is, of course, far from reality while neglecting 
hadron masses, the entropy of the QGP phase can be carried away rapidly by the HHG. 
This looks different for the ideal massless pion gas, where the effective number of hadronic 
degrees of freedom is small, g^ = 3, and consequently, the mixed phase lasts relatively 
long. 
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Our simulation starts as soon as the fireball produced in a central collision of two 
ultra-relativistic heavy nuclei is in local thermal equilibrium. At this initial time tq, we 
have an initial temperature Tq. In the scenario with a phase transition, this temperature 
is assumed sufficiently high for the existence of QGP, which means the hydro dynamical 



description starts in the QGP phase. Inserting the bag model EOS (2.24) in the basic 
differential equation obtained with Bjorken initial conditions ( 2.16| ), one gets the following 
solution for this phase, 

eg(,T) = [e,{To)-B]i^-J +B. (2.38) 

With Eq. ( 2.2lj ), the evolution of temperature follows directly from this solution, 

-^j . (2.39) 

When the decreasing temperature reaches the critical temperature as given by Eq. ( p.33| ), 
the formation of HHG begins. In this mixed phase that lasts until all QGP has hadronized, 



the temperature remains at Tc and the evolution of the energy density is given by Eq. (2.37). 
At the completion of the phase transition, the fireball consists purely of hadronic mat- 



ter. For the pure HHG phase, the basic hydrodynamical equation ( 2.16 ) must be solved 



together with the EOS for the ideal massless pion gas ( 2.29 ) to obtain 

-4/3 



e^{r)=Eh{T^)[^] (2.40) 



and 

n{T) = TA-^\ . (2.41) 

The hydrodynamical treatment is only justified for fluid- like matter. Because of the 
ongoing expansion, the hadronic matter gets more and more dilute. We choose a freeze- 
out temperature Tj as the stopping point of our hydrodynamical simulation. At this 
freeze-out temperature, it is assumed that the thermal interactions in the pion gas vanish 
and that the pions are moving as free particles out of the collision zone. 
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For a closer illustration of the phase transition scenario, we set the model parameters 
to typical values, 

Qq = 37 (two-flavored QGP), 

5/i = 3 (ideal massless pion gas). 

To = 1 fm. 

To = 250 MeV, 

Tc = 170 MeV, 

Tf = 150 MeV. (2.42) 
With these settings, the time scale of the phase transition scenario is determined by the 



entropy density evolution that is shown in Fig. 2.3. Because we assume an adiabatic 



expansion, this graph does not look different for a scenario without a phase transition 



that has the same value of sqTq. In Fig. |2.4| , the proper time dependence of the energy 



density e and the temperature T is illustrated. On the temperature plot, one can see 
clearly the different stages of the collision. As a consequence of the difference in the 
effective number of degrees of freedom, the temperature decrease in the pure QGP phase 
is obviously steeper than in the pure HHG phase. 



Lifetimes of QGP, Mixed and HHG Phase 

In our systematic study of photon yields which will be presented in Chap. ^ we will be 
interested in the contributions from the different collision stages. Therefore, it is important 
to consider also the lifetimes of these stages. The pure QGP phase starts at tq, which is 
a parameter of the model, and ends at 

r? =(!)%« (2.43) 

when the transition temperature Tc is reached according to Eq. ( |2.39| ). The above proper 
time point marks also the beginning of the mixed phase that lasts until A = 0. With 
Eq. ( |2.36| ), one thus gets 

.h _ 9q ( To 



r^ = ^[7^ ro (2.44) 
9h \Tc 
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Figure 2.3: Entropy Density Evolution. The proper time dependence of the entropy 
density is presented for the model parameters, gq = 37, gh = 3, tq = 1 fm, Tq = 250 MeV, 
and Tf = 150 MeV. The values at the onset and offset of the phase transition were obtained 
for Tc = 170 MeV. Ignoring the marks related to the phase transition, this diagram is also 
valid for the no phase transition scenario that has the same value of soTq. Assuming an 
identical thermalization time of tq = 1 fm, the corresponding no phase transition scenario 
must start with the initial entropy density of the above phase transition scenario. This 
demands an initial temperature of To = 578 MeV if the ideal massless pion gas, gh = 3, 
is used to describe the purely hadronic scenario. 
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Figure 2.4: Energy Density and Temperature Evolution in the Phase Transition Scenario. 
The model parameters were set to typical values, gq = 37, gh = 3, tq = 1 fm, To = 
250 MeV, Tc = 170 MeV, and Tf = 150 MeV. 
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for the ending of the mixed phase and the beginning of the pure HHG phase. The pure 
HHG phase ends when the freeze-out temperature Tf is reached which happens at 

3 



9, To , 



(2.45) 



as can be seen from Eq. (2.41). Finahy, using Eqs. ( p. 431) , ( |2.44| ), and ( 2.45| ), one finds for 
the hfetimes the following expressions 



At, 



At/, 



To 



Tc. 



2± 

Qh 



(2.46) 
(2.47) 
(2.48) 



In Chap. 1^, we will refer several times to these expressions and also discuss their tq, Tq, 
Tc, and Tf dependence. 



2.5 The No Phase Transition Scenario 

To decide if a phase transition occurred, one needs to compare experimental data with the- 
oretical predictions. In theoretical calculations, two scenarios always need to be addressed, 
a "phase transition scenario" and a "no phase transition scenario." A clear indication for 
the existence of QGP will be an experimental result that only coincides with theoretical 
predictions obtained in a calculation with a phase transition. The signature cannot be 
considered clean if it can also be deduced in a reasonable calculation without a phase 
transition. Thus, also in this thesis we do not only consider a scenario that exhibits a 
phase transition but also one that is purely hadronic (no phase transition scenario). 

The hydrodynamic description starts at tq with a massless pion gas of temperature Tq 
because no QGP is produced in the no phase transition scenario. Thus, the subsequent 
energy density and temperature evolution is as outlined for the HHG phase of the phase 
transition scenario with the recognition that t^ and need to be replaced by tq and Tq, 
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respectively. To be explicit, we get 



Shir) = e(To) ( - ) (2.49) 



and 

niT)=To(^^y . (2.50) 

As the phase transition scenario, this simulation ends with the freeze-out of the pions at 
temperature Tf. 

For a meaningful direct comparison of a scenario with and another without a phase 
transition, we follow the approach of [^] by assuming identical values of the entropy S 
and the thermalization time tq in both scenarios. Because the entropy is a constant of the 
motion, as is expressed in Eq. ( |2.19| ), this implies identical initial entropy densities s{tq) 
and accordingly also an identical evolution of the entropy density, which can be seen on 
the relation 

where the denominator describes the total fireball volume Vtot(To) as obtained in the 
Bjorken model. Thus, the no phase transition scenario with gf^ effective degrees of freedom 
must have an initial temperature of 

ro^=(^^y^Vo« (2.52) 

to be appropriate for the direct comparison^ with the phase transition scenario that starts 

with a QGP of Qg effective degrees of freedom at Tq. For example, the no phase transition 

scenario that should be compared to the scenario presented in Sec. ^]4| {gg = 37, g^ = 

3, To = 1 fm. To = 250 MeV, = 170 MeV, Tf = 150 MeV) must have the same 

thermalization time tq, and in the case of an ideal massless pion gas, an initial temperature 

of 578 MeV. This extremely high value is a consequence of the small number of effective 

degrees of freedom present in the ideal massless pion gas {gh = 3) and exhibits that the 

ideal pion gas is definitively not the best model for the description of the purely hadronic 

scenario. However, since the emphasis in this thesis is on bremsstrahlung processes in the 

QGP, we now ignore this fact. 

^To be appropriate for the direct comparison with the phase transition scenario means here to start 
with the same initial entropy density as the phase transition scenario. 
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Considering the no phase transition scenario with the following parameters, 



fl'/i = 3 (ideal massless pion gas), 

To = 1 fm. 

To = 578 MeV, 

Tf = 150 MeV, (2.53) 

we find the evolution of entropy density, energy density and temperature as illustrated in 
Figs. |2.3| and |2.5| , respectively. While the entropy density evolution in this purely hadronic 
scenario is identical to the one in the comparable phase transition scenario, the energy 
density and temperature evolution show, of course, a different behavior. 



Lifetime of the HHG Phase 



In the no phase transition scenario, the lifetime of the HHG phase coincides with the 
period in which the fireball is in local thermal equilibrium. This state is reached at the 
initial time tq that marks the starting point of the HHG phase. The ending point is just 
given by the freeze-out temperature Tf, at which the hadrons cease to interact thermally. 
With Eq. (|2]5^), we derive 

(2.54) 



and get for the lifetime of the HHG phase in the no phase transition scenario 

3 



1 



(2.55) 
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Figure 2.5: Energy Density and Temperature Evolution in the No Phase Transition Sce- 
nario. The model parameters were set as follows, gt = 3, tq = 1 fm, Tq = 578 MeV, 
and Tf = 150 MeV. With these settings a direct comparison with the sample scenario 
presented in Sec. ^ (c/g = 37, = 3, tq = I fm, Tq = 250 MeV, = 170 MeV, 
Tf = 150 MeV) becomes meaningful as is explained in the text. 
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Chapter 3 
Photons 



The central element of the research for this thesis were hard thermal photon rates. In 
fact, the advent of new, extraordinary photon rates describing bremsstrahlung processes 
in QGP triggered the investigation presented in this thesis. This chapter will review the 
present-day rates for hard thermal photons from QGP and HHG. 

As mentioned in the introduction, photons are interesting probes of the fireball. Be- 
cause of their purely electromagnetic nature, they have mean free paths much larger than 
the transverse size of the fireball, Xmfp ^ o". Thus, most photons produced in the fireball 
do not suffer any final state interaction. They reach the experimental detector directly 
from their origin in space-time. As a consequence, photon distributions stay far from 



equilibrium |21| and provide undisturbed footprints from the various stages of nuclear 
matter in which they have been produced. Photons seem in this sense promising for the 
confirmation of QGP possibly produced in ultra-relativistic heavy ion collisions. 

High-energy nucleus-nucleus collisions exhibit many (potential) photon sources. There 
are hard scattering processes of the initial partons during the very early stages of the 
collision. Among these processes that bring the system locally into thermal equilibrium 
are also photon producing reactions. Photons from these processes are referred to as 
prompt photons. When local thermal equilibrium is reached, there will be thermal photon 
production. Because the electrically charged constituents that emit photons are different 
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in QGP and HHG, it is sensible to consider thermal photon emission from QGP and HHG 
separately. After freeze-out, hadronic decay processes as the Dalitz decays, — > 77 and 
77 — > 77, contribute significantly in the measured photon spectra. 

In this thesis, we concentrate on thermal photon production, so prompt and decay pho- 
tons are considered background. Because experiments provide direct photon yields, decay 
contributions are already subtracted as will be explained in Chap. |5[ The situation looks 
different for prompt photons that are part of the direct photon spectrum. It is the task of 
theorists to calculate and subtract the prompt photon yields in order to get the thermal 
photon spectrum. However, since the photon production in pre-equilibrium scatterings 
is mainly important for the very high momentum region of the considered spectra [2 4( 1, 
we only discuss prompt photon production qualitatively. Without the quantitative treat- 
ment, one must keep in mind that there are photon emitting processes besides thermal 
photon production when we inspect the direct photon yields at very high momenta. 

This chapter is organized as follows. We start with a brief discussion of prompt 
photons. Next, a detailed review of the most recent thermal photons rates for QGP and 
HHG is presented. Decay photons are not covered because they do not appear in the 
experimental direct photon data. A description of how we get the photon spectrum from 
the production rates concludes this part of the thesis. 



3.1 Prompt Photons 



In the very early pre-equilibrium stage of an ultra-relativistic heavy reaction, many hard 
scattering processes between the partons of the initial colliding nuclei steer the system 
towards local thermal equilibrium. Among these hard processes, there are also photon 
producing reactions, such as Compton scattering {qg qj, qg — > ^7), gg-annihilation 
{qq — > 57) and higher-order bremsstrahlung processes {qg qgj, qg — > qg'j, etc.). 
Because the momentum transfer in these processes is large, the elementary rates on 
the parton level, {dN / d? p^) ab^-yc-, can be calculated reliably within perturbative (zero- 
temperature) QCD. These rates must then be folded with the structure functions F of 
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the colliding nuclei A and B, 

as is derived in [^]. In the above expression, C stands for the final state particle(s) 
on the hadronic level produced besides the photon, c denotes the same quantity on the 
partonic level, and Fjj{xj)dxj is the probability of having a parton j inside the nucleus 
J, where the parton carries a momentum fraction between xj and Xj + dxj of the nucleus 
momentum. By summing the rates of all processes AB 'jC up to some specified 
order in the coupling constants incoherently, one obtains the total prompt photon yield 
dN / Sp^ in the specified order. 

The uncertainties in the extraction of the prompt photon yields are dominated by two 
aspects connected to the nuclear structure functions Fj j{xj). First, even for relatively 
hard photons (-E^ ~ 4 GeV), one is in the low-x region, where the form of the structure 



functions is a matter of current research 3£]. Second, the structure functions for 
the nuclei cannot be simply extrapolated from the nucleon structure functions: one has 
to take into account medium effects that are also the subject of ongoing discussions. 
Here, we do not go beyond the above qualitative remarks on prompt photons. For a 
quantitative treatment, the reader is referred to a recent investigation with emphasis on 
nuclear shadowing effects [40|. 



3.2 Thermal Photons 

The thermally moving constituents of the strongly interacting continuum, produced in an 
ultra-relativistic heavy ion reaction, undergo reactions and, in the case of hadrons, also 
decays. These processes are governed by the thermal distributions of the participating 
particles. Because many constituents of the continuum carry electrical charge, many 
reactions produce photons. These are the photons we refer to as thermal photons. 
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3.2.1 Thermal Photons from Quark-Gluon Plasma 



Thermal photon production from QGP is examined in finite temperature QCD [41|. The 
lowest order photon emitting reactions in the QGP are gg-annihilation, 

qq gj, (3.2) 

and Compton scattering with an initial gluon, 

qg qj, (3.3) 

qg 97, (3.4) 



which are shown as Feynman diagrams in Fig. 3.1. By considering these processes for 
massless quarks in QCD, a logarithmically divergent behavior of the production rates is 
encountered for soft momentum transfers. This infrared divergence is well known since it 
appears in every QED and QCD tree diagram process in which a soft, massless particle is 
exchanged. However, QGP medium effects result in a non-vanishing thermal quark mass 
that serves as an infrared cutoff and renders the production rate being finite. A systematic 
treatment of the QGP medium effects that does not spoil gauge invariance became avail- 
able with the Braaten-Pisarski method known as hard thermal loop (HTL) resummation 
technique [^] . Using this technique, the production rates of hard {E^ S> T) thermal pho- 
tons from the above lowest order reactions have been calculated for thermal and chemical 
equilibrium ^] and also for chemical non-equilibrium [^9| , where deviations from the 
thermal distributions were described by fugacities. A close inspection of photon spectra 
from a fireball in chemical non-equilibrium has been found as being important but it is 
beyond the scope of this study. We concentrate only on photon emission from a system 
in thermal and chemical equilibrium at zero baryon density. The case of finite baryon 



density, which has also been investigated |51], is here of no relevance. 



How does one proceed (more explicitly) in computing the production rate of hard 
thermal photons from gg-annihilation (qq — > g^) and Compton scattering with an initial 



gluon (qg — > qj, qg — > ^7)? Following the Braaten-Yuan prescription |43], a parameter A 
is introduced that separates high momentum transfers of the order T from low momentum 
transfers of the order gT, where the strong coupling constant is assumed much smaller 
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Figure 3.1: gg- Annihilation and Compton Scattering with an Initial Gluon. The Feynman 
diagrams display the lowest order reactions in which a real photon is produced: qq- 
annihilation (left hand side) and Compton scattering with an initial gluon (right hand 
side). Only the Compton scattering diagrams for a quark are shown because the same 
topology with reverted particle flow holds for the Compton scattering of an anti-quark. 
For massless quarks and soft momentum transfers, each process exhibits the well known 
logarithmically divergent behavior. In QGP, however, medium effects provide the quarks 
with a finite thermal mass that regulates all infrared divergences. A consistent inclusion of 
the medium effects is achieved within the Braaten-Pisarski method, where the soft quark 
propagator is replaced by an effective propagator that emerges from the resummation of 
hard thermal loops. 
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than one, <C 1. This decomposes the rate into a hard and a soft part. The hard part has 
the infrared cutoff A and is calculated in relativistic kinetic theory with bare propagators 
and bare vertices which means that no medium effects are taken into account. The thermal 
distributions enter the derivation of the hard part in the following way 

X ^ |A^|2 {27r)^6Hpi +P2-P3- P-y), (3.5) 

where l-^P is the squared scattering amplitude of the considered process summed over 
the initial and final parton states. The particles in the initial state are labeled by 1 and 
2, while those in final state are labeled by 3 and 7. The /'s are the Fermi-Dirac or 
Bose-Einstein distribution functions depending on the spin of the corresponding particle. 
A fermion in the final state is Pauli suppressed while a boson in the final state is Bose 
enhanced. For the soft part with the ultraviolet cutoff A, medium effects are crucial 
and any soft bare quark propagator is replaced by the HTL resummed propagator that 
contains the finite thermal quark mass. The soft quark propagators in other words become 
dressed. The soft part of the thermal emission rate is calculated from the retarded photon 
self-energy according to the relation ^ 

where n^''^ is the retarded polarization tensor of the photon and /b is the Bose-Einstein 
distribution function, (/^-annihilation and Compton scattering are contained in the one- 
loop contribution to the photon self-energy with one quark propagator dressed]^. This can 



be seen by applying the thermal cutting rules 47] on the diagram shown in Fig. 3.2 . 
By summing the hard and the soft part, both of which depend on A, a finite result is 
obtained, which is independent of the separation-parameter A. 



^The one-loop contribution to the photon self-energy with both quark propagators dressed does not 
enter into the computation of the hard thermal photon rates because the photon cannot be hard if both 
quarks in the loop are soft. 
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Figure 3.2: Effective One-Loop Contribution to the Plioton Self-Energy. The hard thermal 



photon production rates for the processes shown in Fig. |3.1| arise from the imaginary part 
of this diagram by applying thermal cutting rules. One quark in the loop can be soft 
which means that its bare propagator must be dressed to avoid an infrared singularity. 
The dressing is indicated in the black blob on the lower quark propagator. 

The net rate in the limit of hard photons, Ey ^ T, has the following form 
dN 



E 



l-loop f ^ 



where c = 0.23 is a constant, is the electric charge of the quark with flavor / in units 
of the electron charge e, and the sum runs over all flavors assumed in the QGP, e.g., 

for a two-flavored QGP. We use expression (^3) with the standard value A'^c = 3 and 
Cp = 4/3, which is the corresponding Casimir operator of the fundamental representation 
of color SU(3). Further, we insert for the QED coupling constant a = 1/137 and for the 



QCD coupling constant the lattice QCD data parameterization |52] 



6 vr 

"^^^^ " (33-2iVj) H8T/T,)- ^^-^^ 



The rate (3.7) is labeled with the subscript 1 — loop to remind that the soft part was 
computed from the one-loop contribution to the photon self-energy. Its dependence on 
the photon energy E-f for two fixed temperatures, T = 150 MeV and T = 250 MeV, is 



presented with the dashed lines in Fig. 3.3, where a two-flavored QGP has been assumed. 



The 1 — loop rate has been employed in many investigations on real photons as the only 
measure for thermal photon production in QGP |23, 24, 26, 37, 4C]. This must he modified 
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Figure 3.3: Thermal Production Rates of Hard Photons from QGP and HHG. The 

Ey dependence of the thermal photon production is presented for two temperatures, 
T = 150 MeV and T = 250 MeV, where the dashed, dot-dashed, and solid lines rep- 
resent the contributions from 1 — loop, bremss, and qq — aws processes in a two-flavored 
QGP , respectively, and the dotted line represents the contribution from np — and 
Tvp ^ ai ^ 7r7 processes in HHG. Because all rates have been computed in the limit of 
hard photons, Ej S> T, the shape of the curves for Ej < 1 GeV should be ignored. It is 
important to notice that already the static thermal photon spectra exhibit the qq — aws 
process as the dominating one for E^ > 1 GeV. This is a remarkable result since the 
qq — aws process is of higher order than the 1 — loop processes. 
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Figure 3.4: Effective Two-Loop Contribution to the Plioton Self-Energy. In thermal 
QCD, bremsstrahlung processes have been investigated on the two-loop contributions to 
the photon self-energy within the framework of the Braaten-Pisarski method. By applying 
thermal cutting rules on the diagrams, the bremsstrahlung processes that are illustrated in 
Fig. emerge. The black blob indicates an effective gluon propagator which is necessary 
since the gluon can be soft. It is basically this effective gluon propagator which causes 
the fact that bremsstrahlung processes arise in the same order of the coupling constants 
as the one-loop processes. 

because bremsstrahlung processes considered in thermal QCD turned out to contribute at 
the same order in the coupling constants as the one- loop processes [^0| . The thermal rates 
for the bremsstrahlung processes have been calculated from the two-loop contributions to 



the photon self-energy that are pictured in Fig. |3.4| . The thick black dot on the gluon 
propagator indicates an HTL resummed propagator which is necessary since the gluon 
can be soft. By applying the thermal cutting rules on the two-loop contribution to the 
photon self-energy, physical scattering processes such as the ones shown in Fig. are 
obtained. It can be seen that besides "ordinary" bremsstrahlung. Fig. |3.5| (a), another 
process of gg-annihilation with an additional scattering in the medium. Fig. |3.5| (b), is 
contained in the two-loop contributions to the photon self-energy. This process has in fact 
been found as the dominating source for hard photons in the QGP. 
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Figure 3.5: Feynman Diagrams of the Bremsstrahlung Processes. By applying thermal 
cutting rules on the effective two-loop contributions to the photon-self energy, one finds 
the physical bremsstrahlung processes with a real photon in the final state. The diagram 
on the left hand side is just the ordinary bremsstrahlung process involving two quarks. 
The dressing on the gluon propagator indicated by the black blob is necessary to avoid an 
quadratical infrared divergence for soft momentum transfers. There are similar processes 
which are not shown but also described in the production rate bremss: qq — »• qq^y, qq — »• 
QQl^ Q9 ~^ Q97j qg — ^ qgj. The diagram on the right hand side is a reaction that can 
be described as gg-annihilation with an additional scattering on a quark. Again, other 
similar processes obtained from the two-loop contributions to the photon-self energy are 
not illustrated but included in the production rate qq — aws, such as qqq 57 and 
qqg — > 57. Reactions of the type shown on the right hand side, qq — aws, are the ones 
that surprisingly dominate the thermal photon spectra from the QGP. 
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The technical derivation of Aurenche et al. leads in the limit E"^ ^ T to the production 
rate for ordinary bremsstrahlung photons 



E 



7 



= 15 ln(2)( Jt - J^)^ ( ^ e^) aa, e"^-/^ (3.10) 

bremss f 



and to the production rate for photons from gg-annihilation with an additional scattering 
in the medium 



= ^^J---'-)^{T. 4) T (3-11) 

qq—aws f 



where the factor {Jt — Jl) in the above expressions stands for the difference of two integrals 
that depend only on the number of colors N^. and flavors Nf assumed in the QGP. With 
the standard value Nc = 3, one gets for a two-flavored QGP, 

{Jt - Jl) = 8.71. (3.12) 

The dependence of ( |3.10| ) and ( |3.11| ) on the photon energy E^ for two fixed temperatures, 
T = 150 MeV and T = 250 MeV, under the assumption of a two-flavored QGP is shown 
graphically in Fig. with the dot-dashed lines representing the bremss contribution 
and the solid lines representing the qq — aws contribution. The dominance of the qq- 
annihilation process with an additional scattering in the medium can be seen clearly. It 
can be traced back to the factor E^T exp{—E^/T) that favors the production of hard 
thermal photons. In Chap. ^, one will realize also on the thermal photon spectra that the 
bremsstrahlung processes, bremss and qq — aws, do not only give small corrections but 
lead to a significant enhancement in the thermal photon yield from QGP. 

There are several uncertainties in the thermal photon production rates which have 
been derived in finite temperature QCD by applying the Braaten-Pisarski method to 
account for long range medium effects in the QGP. The rates listed explicitly have been 
obtained for thermal and chemical equilibrium. As is shown in |48, 49, pO], abandoning the 



assumption of a chemically equilibrated QGP significantly alters the one-loop result. Thus, 
an inspection of bremsstrahlung processes from a system in chemical non-equilibrium 
would be very interesting and, in fact, it should be the subject of a future systematic 
investigation. Another uncertainty is connected to processes contained in the three or 
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more loop contributions to the photon self-energy. Photon self-energy terms with a higher 
number of loops could again contribute in the same order of the coupling constants as 
the one- and two-loop contributions. A confirmation of this speculation would definitively 
raise fundamental questions concerning the validity of the effective perturbative expansion 
based on the HTL resummation technique. Finally, the HTL resummation technique 
is based on the assumption g <^ 1 which does not hold even in the hot stage of an 
ultra-relativistic heavy ion collision: for a two- flavored QGP at a high temperature, e.g., 
T = 3Tc, Eq. ( |3.9[ ) gives Os = 0.2 which corresponds to g = -v/ivraj = 1.6. However, an 
effective field theory cleaner than the one achieved from the HTL resummation technique 
and also suitable to calculate thermal photon rates for QGP is not available. 

3.2.2 Thermal Photons from Hot Hadronic Gas 

In our simple dynamical model of nuclear collisions, we treat the HHG as an ideal massless 
pion gas. This is, of course, only an approximation in which thermodynamic quantities are 
obtained analytically. Pions are not massless and the HHG phase will also exhibit other 
mesons, such as p, r], and to mesons. We do not account for these facts in the consideration 
of the fireball evolution as already explained in Chap. |2|. However, we implement thermal 
photon production rates from reactions of mesons more massive than pions. This approach 
is to some extend inconsistent, but sufficiently precise for the purpose of this investigation 
that centers on the effects of the bremsstrahlung processes in QGP. 

The thermal photon rates for reactions of the thermalized secondary hadrons have 
been derived in an effective field theory that is renormalizable and provides gauge invariant 
results. Scattering processes of vr, p, and r] mesons, 



TTp 7r7, (3.13) 

TTTT /97, (3-14) 

TTTT — > ?77, (3.15) 

TTT] 7r7, (3.16) 

TTTT 77, (3-17) 
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and also vector- meson decays, 



p° ^ 7r+7r-7, (3.18) 
uj 7r°7, (3.19) 

have been examined and the Compton scattering of a pion on an initial p meson, vrp n^, 



has been explored as the dominating photon source in HHG |21|. A later investigation of 



the np — > 7r7 scattering but this time through the ai-resonance, 

vrp — > ai — > TTj, (3.20) 
discovered a significant enhancement of the np contribution |27]. Concentrating 



on hard thermal photons, we consider the dominating processes (|1|) and (|3^ ) as the 



only photon sources in the HHG and apply the parametrization from [p7|] divided by a 



factor of two, which was suggested by Kevin Haglin |28| due to an encountered isospin 
miscounting in [^], 



E-y 



dN 



had 



where E^ and T should be given in units of GeV. The dotted lines in Fig. 3^ illustrate the 
dependence of this rate on the photon energy E^ for two fixed temperatures, T = 150 MeV 
and T = 250 MeV. Other decay processes besides the ones listed above can be neglected: 
either they are too slow, such as the Dalitz decays, or they involve more massive and 
consequently stronger Boltzmann-suppressed hadrons. 

Uncertainties arise mainly from the assumption of a thermally and chemically equi- 
librated HHG phase and from potential hadronic in-medium effects, such as decreasing 
masses and increasing widths, that have been ignored in the derivation of the above pa- 
rameterization. In addition, there is no fundamental theory for the strong force on the 
hadron level which means that the coupling constants of the effective theory must be 
inferred from various experiments that do not display the environment produced in an 
ultra-relativistic heavy ion collision. 
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3.3 Photon Spectra 



It has been mentioned before that photon rates must be convoluted with the space-time 
evolution of the fireball to obtain the photon spectra comparable with experimental data. 
This section details on the explicit form of the integration ( |2.2D in the fireball picture 
described. 

Because we concentrate on central collisions with vanishing impact parameter and 
consider an only longitudinally expanding tube of strongly interacting matter, the inte- 
grations in the directions transverse to the beam can be performed independently of the 
production rate 

' d^x = TTR\ f dtdz, (3.22) 



where expression ( p. 11 ) is used. In the remaining two-dimensional integration, we sub- 



stitute the center-of-mass coordinates t and z of the emitting fluid cell by proper time r 
and rapidity fluid cell rapidity y', 

t = rsinhy', (3.23) 
z = r cosh 2/', (3.24) 

since in the latter variables one takes advantage of the symmetries contained in the Bjorken 
model. The fluid cell rapidity has been renamed from y to y' and in the following y will 
be referred to as the photon rapidity in the center-of-mass frame. The substitution results 
in the form 

rT2 r + Vnucl 

dtdz= dTT dy', (3.25) 

where ri and T2 denote the starting and ending point of the considered collision phase, 
respectively. If any nuclear stopping is neglected, ynud is the center-of-mass rapidity of 
the projectiles and the identity 

Vnud = arcosh {^ 2-A^Ge\ ) ^^''^^^ 
holds, which means ^/s limits the range of the fluid cell rapidities. Next, we rewrite 

d^p d'^p_j_ dy'' ^ ^ 
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where E, p, pj_, and y represent energy, momentum, transverse momentum, and rapidity 
of the thermal photons, respectively. The index 7, formerly used to label a quantity that 
refers to the emitted thermal photon, has been suppressed. By inserting Eqs. ( |3.22 ), 



( 3.25| ), and (3.27) into the integration (p.2|), one arrives finally at the expression that is 



evaluated numerically to obtain the thermal photon spectra in the center-of-mass (CM) 
frame 

,9 r =t^Ra dTT dy' [E ,„ . (3.28 

J^P±dyJ(^CM) ^ri J-Vnuci V d'^xd'ipj ^LR) 

The integrand of the above two-dimensional integral is the production rate in the consid- 
ered collision phase (pure QGP, mixed, or pure HHG phase). Because production rate 
expressions have been computed for the local rest frame (LR) of the emitting fluid cell, 
the following form of photon energy must be used in the integrand, 

%R) =P± cosh(y'-y). (3.29) 

During the mixed phase, the QGP volume fraction A enters the integrand, 

dN \ / dN \^^^ r .1 / dN \"^^ 

dHd-^pJi^LR) V d^xd^pJ^LR) V d'^x d-^p J i^LR) 

with the rates labeled by the superscripts QGP and HHG referring to the QGP and the 
HHG state of matter, respectively. 
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Chapter 4 



Systematic Investigation 



After the introduction on ultra-relativistic heavy ion physics, the description of the simple 
model for heavy ion reactions in the transparent region, and the discussion of the thermal 
photon production rates, we now present the results of the systematic investigation on the 
thermal photon spectra. The starting point are the spectra obtained in the two sample 
scenarios that have been introduced in Chap. ^, one with phase transition and the other 
without. Then, we concentrate only on the phase transition scenario, where the one 
presented in Chap. ^ serves as the basis for the systematic investigation, in which the 
dependence of the spectra is carefully examined on the mass number of the projectiles A, 
the projectile rapidity VnucU the thermalization time tq, the initial temperature Tq, the 
transition temperature Tc, and the freeze- out temperature Tf. 



4.1 Thermal Photon Spectra in the Sample Scenarios 



The fireball evolution provided by the simple, well understood model explicitly described 
in Chap. ^ was illustrated on two sample scenarios, a phase transition scenario and a no 
phase transition scenario. The model parameters were specified for the phase transition 
sample scenario in ( 2.42| ) and for the no phase transition sample scenario in (2.53), where 
Tq was adjusted according to Eq. ( p. 52 ) to allow a direct comparison of both scenarios. We 
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now extract the thermal photon spectra for both sample scenarios by applying Eq. ( p. 28] ) 
with the necessary additional parameters set to 



Vnud = 8.6, (4.1) 

corresponding to the pursued LHC center-of-mass collision energy of y/s = 5500 A ■ GeV, 
and 

A = 208, (4.2) 

corresponding to the mass number of the lead nuclei that will be the projectiles in the 
LHC heavy ion physics program. To study photons from the central, low-baryon density 
region which is best described by the employed model, all photon spectra are presented 
for photons at mid-rapidity, 

y = 0. (4.3) 



4.1.1 Thermal Photon Spectra in the Phase Transition Scenario 

In the phase transition scenario, several combinations of the spectra give insights into the 
importance of the different rates, the two states of matter (QGP and HHG), and the three 
collision phases (pure QGP phase, MP, and pure HHG phase), where we are primarily 
interested in the effects of the QGP bremsstrahlung processes. The upper diagram in 
Fig. |4.1| shows therefore the totalQ thermal photon yields with bremsstrahlung processes 
(solid line) referred to as New and without bremsstrahlung processes (dashed line) referred 
to as Old. The effect of the QGP bremsstrahlung processes is clearly most pronounced for 
high-p_|_ photons. This behavior is due to the fact that the high-p_|_ region is mainly popu- 
lated through processes at high temperatures. In the phase transition scenario, these are 
the processes in the pure QGP phase and the MP to which bremsstrahlung processes (if 



included) belong. The lower diagram in Fig. iA illustrates the decomposed contribution 
of these processes from the QGP state of matter, where the dashed, dot-short-dashed, and 
dot-long-dashed lines indicate the contributions from the 1 — loop, bremss, and qq — aws 
production rates, respectively, and the solid line indicates the sum of the contributions 



^total — pure QGP phase contribution + MP contribution + pure HHG phase contribution 
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from the three QGP rates. This diagram exhibits most clearly the astonishing importance 
of the bremsstrahlung processes: the inclusion of the bremsstrahlung processes enhances 
the total thermal photon yield from the QGP state of matter by about one order of mag- 
nitude over the complete considered pj_-range, and already the qq — aws processes alone 
outshine the 1 — loop processes that give the total yield from the QGP if bremsstrahlung 
processes are neglected. 

In Fig. the yields from the two states of matter, QGP and HHG, can be compared 
with (New) and without (Old) the QGP two-loop contributions taken into account. The 
spectra from QGP and HHG are presented in the dashed and dotted lines, respectively, 
and the sum is presented in the solid line. Without bremsstrahlung processes in the QGP, 
HHG controls the total thermal photon spectrum almost over the complete considered 
transverse momentum range, 1 GeV < p± < 4.5 GeV, while with these processes, QGP 
becomes the dominant photon source for p± > 3 GeV. 

It is also interesting to inspect the photon yields from the different stages of the 
collision, the pure QGP phase, the MP, and the pure HHG phase, that are shown in 
Fig. The upper plot, where the bremss and the qq — aws rates were included, exhibits 
the pure QGP phase (dashed line) as the dominant thermal photon source in the high-p_L 
region. Further, in the upper and the lower plot, where for the QGP state of matter only 
the 1 — loop rate was implemented, one can clearly identify different slopes corresponding 
to different temperature ranges of the specific collision stage. Due to the Boltzmann factor, 
exp(—E/T), that is contained in each of the applied rates as the essential factor for the 
shape of the spectrum, the relation between slope and temperature can be expressed 
approximately in the form 




1 



(4.4) 



which means a flat spectrum indicates a high temperature. 
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Figure 4.1: Total Thermal Photon Spectra and Thermal Photon Spectra from the QGP 
State of Matter in the Phase Transition Sample Scenario. The total thermal photon 
yields are shown in the upper plot, where the solid and dashed lines are obtained with 
and without QGP bremsstrahlung processes, respectively. In the lower plot, the thermal 
photon yields from the QGP state of matter are presented, where the dashed, dot-short- 
dashed, and dot-long-dashed lines indicate the spectra from the 1 — loop, bremss, and 
qq — aws rates, respectively. The total QGP contribution is given in the solid line if 
bremsstrahlung processes are included. 54 





Figure 4.2: Thermal Photon Spectra from the QGP and the HHG State of Matter in 
the Phase Transition Sample Scenario. The thermal photon yields from the QGP and 
the HHG state of matter arc represented in the dashed and the dotted lines, respectively, 
while the sum of both is the total thermal photon yield represented in the solid line. QGP 
bremsstrahlung processes are included in the upper plot and neglected in the lower plot. 



55 



Photon Spectrum (QGP - MP - HHG) New 



= total _ 

totalg 




1 1.5 2 2.5 3 3.5 4 4.5 5 



P± [GeV] 



Photon Spectrum (QGP - MP - HHG) Old 




1 1.5 2 2.5 3 3.5 4 4.5 5 

P± [GeV] 

Figure 4.3: Thermal Photon Spectra from the Different CoUision Stages in the Phase 
Transition Sample Scenario. The dashed, dot-dashed, and dotted lines correspond to the 
thermal photon yields from the pure QGP, the MP, and the pure HHG stages of the 
collision, respectively, and the solid line corresponds to the total thermal photon yield, 
where QGP two-loop processes are included in the upper diagram and neglected in the 
lower diagram. 
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Photon Spectrum (HHG) 
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Figure 4.4: Thermal Photon Spectrum in the No Phase Transition Sample Scenario. 



4.1.2 Thermal Photon Spectra in the No Phase Transition Scenario 



While in the phase transition scenario it was interesting to analyze the total photon spec- 
trum in several ways, the basic message from the no phase transition scenario can be 



illustrated in just one line which is presented in Fig. 4.4. This line, which indicates the 



thermal photon yield from the HHG state of matter, is obtained with the production 



rate given in Eq. (3.21). It is trivial that the scenario without a phase transition, in 
which no QGP state of matter is assumed, does not allow to study the effects of QGP 
bremsstrahlung processes. Nevertheless, the comparison of the total thermal photon spec- 
trum from the no phase transition scenario with the one from the phase transition scenario 
is very important because the quality of thermal photons as a potential signature for the 
QGP formation can be tested in this way. Here, the comparison of the total thermal 
photon spectra exhibits a significant abundance in the no phase transition scenario that 
increases from one order of magnitude for photons with p± = 1 GeV up to about six or- 
ders of magnitude for photons with p± = 5 GeV. This behavior can be traced back to the 



temperature evolution shown in the lower half of Fig. 2.5. The high initial temperature 



of To = 578 MeV in the no phase transition scenario, which was chosen to have the same 
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initial entropy density as in the phase transition scenario, leads to a very high mean tem- 
perature that benefits thermal photon production especially in the high-px range. In fact, 
the extremely high value for Tq was a consequence of the massless pion gas used to model 
the HHG state of matter and its small number of effective degrees of freedom (g^ = 3) that 
entered the derivation of the initial temperature from the initial entropy density. A higher 
number of effective degrees of freedom, which is present in a more realistic description 



of the HHG including finite hadron masses and a higher number of hadrons |26|, would 
result in a lower initial temperature and consequently in a lower thermal photon yield. 
For an adequate treatment of the purely hadronic scenario, one must therefore implement 
such a more realistic description of the HHG that includes finite hadron masses and a 
higher number of hadrons. This will not be presented here and the important careful 
comparison of the thermal photon spectra from the phase transition scenario with those 
from the no phase transition scenario is postponed for future work, which means that we 
will concentrate from now on only on the scenario with a phase transition. 

The remaining part of this chapter is devoted to the dependence of the thermal photon 
spectra on the model parameters in the phase transition scenario. We take the phase 
transition sample scenario with the parameter settings (2.42), ( f4.1[ ), and ( f4.2[ ) as the 



basis. Then, sticking to the two-flavored ideal massless parton gas {gq = 37) and to 
the ideal massless pion gas {gh = 3), we vary every other parameter separately over an 
instructive range and inspect the corresponding influence on the thermal photon spectra. 



4.2 A - Mass Number of Projectile 

The mass number A of the projectile governs the transverse size of the fireball tube. Since 
we concentrate on central collisions with zero impact parameter, the transverse size of the 
only longitudinally expanding fireball equals the transverse size of the projectile. The 
thermal photon spectra are therefore proportional to the geometrical cross section of the 
projectile, 

'^^ oc 7ri?i, (4.5) 



d?p_\_ dy 
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Figure 4.5: Total Thermal Photon Spectra for A = 32 (Upper Curves) and A = 197 



(Lower Curves) in the Phase Transition Scenario. As in the upper diagram of Fig. 4.1, 
the solid lines are obtained if QGP bremsstrahlung processes are included and the dashed 
lines are obtained if QGP bremsstrahlung processes are neglected. 



which can be seen from Eq. ( ^3.28 ). This proportionality together with the phenomeno- 
logical formula (2.11) for the nuclear radius leads directly to the ^-dependence of the 
thermal photon spectra, 

'^^ oc A'/'. (4.6) 



(Pp± dy 

In Fig. the total photon spectrum is presented for sulfur {A = 32) and gold {A = 197) 
projectiles with (solid line) and without (dashed line) bremsstrahlung processes included. 
This diagram illustrates clearly the dependence of the thermal photon spectra on the mass 
number A of the projectile. With the above expression ( [4.6D , it is found that the displayed 
spectra differ by a factor of about 3.5. For lead {A = 208) projectiles, the total thermal 
photon spectra have already been shown in the upper half of Fig. and by comparing 
these spectra with the ones for the gold projectiles, only a marginal difference is found. 

Different accelerators provide different projectiles depending on the employed heavy 



ion injector. Table 4.1 summarizes the projectiles that have been applied at the CERN 
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Accelerator / 


Beam- Target / 


Experiment 


A 




Vnucl 


dN/dy 


To 


Collider 


Beam-Beam 






[A-GeV] 






[MeV] 


SPS 




WA80 


32 


20 


3.0 


200 


185 




208 _208p^ 


WA98 


208 


17 


2.8 


800 


195 


RHIC 


197^^ _197^^ 


PHENIX 


197 


200 


5.3 


>1200 


>225 


LHC 


208 _208p^ 


ALICE 


208 


5500 


8.6 


>2500 


>280 



Table 4.1: Accelerator- and Experiment Specific Quantities Relevant for the Extraction of 
the Thermal Photon Spectra. The entries for the projectile rapidity Unud were obtained 
with Eq. (|3.26D from the corresponding center-of-mass energy ^/s. For the multiplicity 
distribution dN/dy, we specified values in accordance with the literature [|^, |5^, ^ [5^ . 



With these values, the Tq entries were computed using Eq. (4^) and the canonical value 
for the thermalization time tq = 1 fm. 



SPS and that will be applied at the BNL RHIC and the CERN LHC. Collisions of gold 
and lead nuclei support the highest thermal photon yields but unfortunately also the 
highest background yields. However, due to the progressive photon spectrometer LEDA 
used in the WA98 experiment, the direct photon data from WA98 is expected to provide 
interesting information soon on the thermal photon production in ultra-relativistic heavy 
ion reactions. 



4.3 ynuci - Projectile Rapidity 

The projectile rapidity ynud limits the rapidity range of the photon emitting fluid cells 
if any stopping of the projectiles during the collision is neglected. It is calculated from 
the center-of-mass energy according to Eq. ( ^.26| ), which has also been used for the com- 



putation of the accelerator specific values of ynud listed in Tab. 4.1. The assumption 
of no nuclear stopping in an ultra-relativistic heavy ion collision is, of course, arguable. 
However, the dependence of the thermal photon spectra at mid-rapidity on the fiuid cell 



rapidity limits, or ynud-, is anyway extremely weak. Figure i.t illustrates this fact on the 



emission of mid-rapidity photons, ?/ = 0, with two transverse momenta, p_L = 1 GeV and 
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Photon Radiation at T = 250 MeV 
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Figure 4.6: Thermal Emission of Mid-Rapidity Photons with p_L = 1 GeVandpj_ = 5 GeV 
as a Function of the Fluid Cell Rapidity y' at the Temperature T = 250 MeV. As in 



Fig. 3.3, the dashed, dot-dashed, solid, and dotted lines indicate the thermal photon 
production rates 1 — loop, bremss, qq — aws, and had, respectively. 

p± = 5 GeV, from fluid cells with different rapidities y' at the temperature T = 250 MeV. 
It can be seen that the contribution to the mid-rapidity photon yield decreases rapidly 
for both transverse momenta with increasing fluid cell rapidity y' . This behavior is easily 
understood as a consequence of the Boltzmann factor, which is an essential part of every 
thermal photon production rate considered. For example, a fluid cell at the rapidity y' = 4 
must emit a photon having an energy of ^(l_r) = 1 GeV • cosh 4 ~ 27 GeV for it to appear 
at mid-rapidity with a transverse momentum of p± = 1 GeV. Due to the extremely high 
energy, the Boltzmann factor makes this event an extremely rare one. In conclusion, we 
can neglect the weak dependence on ynud ^ we are concentrating on the photon spectra 
at mid-rapidity. 
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4.4 Initial Conditions 



The initial conditions determine the complete evolution of the strongly interacting con- 
tinuum and, therefore, their significant influence on the thermal photon spectra is clear in 
advance. Properties of the collision, as the mass number A of the projectile, the projectile 
rapidity i/nuch the impact parameter b, should unambiguously specify these important 
initial conditions, but there are uncertainties as the nuclear stopping power and the time 
scale for equilibration that complicate their assessment substantially. How can we then 
estimate the initial conditions? A first hint on the thermalization time is given by the 
strong interaction (proper) time scale [p^] , 

TO ^ ^ 1 fm, (4.7) 

^QCD 

which coincides with the canonical value, tq = 1 fm, used in many investigations. In the 
Bjorken model, an additional expression can be derived that relates the thermalization 
time Tq and the initial temperature Tq to the multiplicity distribution dN/dy [^] 

3 _ 27r^ dN 

^°^°-45C(3)vri?i4a, rf^' ^^'^^ 

where for the phase transition scenario with a two-fiavored QGP in the initial state, 
= cLq = 377r^/90. While the multiplicity distribution dN/dy at the CERN SPS can be 
measured experimentally, it must be extrapolated from high energy p-p and p-p collision 
data for RHIC and LHC. By going through the literature [|^, we could specify 



the dN/dy values listed in Tab. iA, where the RHIC and LHC entries present lower 



limits. The initial temperatures given in this table were calculated with Eq. ([4.8[) for the 



canonical thermalization time, tq = 1 fm. In several works |24, 37 1, Eq. ( [4. 8]) has been 
used together with another relation for tq and Tq [^] 

ro^:^ (4.9) 
61q 

and with measurements or estimations of dN/dy as the key to the initial conditions. A 



different approach in the framework of hydrodynamics is presented in |25], where hadronic 
spectra are fitted to experimental data by fixing the initial conditions, which are then used 
to extract photon and dilepton yields. More sophisticated theoretical calculations beyond 
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hydrodynamics are also available to determine the initial conditions since transport models 
can simulate the early non-equilibrium phase of the reaction from kinetic theory. In 
transport models on the partonic level, e.g., the parton cascade model (PCM) parton 
momentum distributions are considered locally and the corresponding cell is revealed as 
being in local thermal equilibrium if these distributions exhibit an exponential (thermal) 
and isotropic distribution. Here, no initial conditions will be derived from observables or 
transport models, we instead restrict ourselves to the variation of tq and Tq within the 
realistic ranges that contain values from the methods discussed above. 



4.4.1 To - Thermalization Time 

The thermalization time tq indicates the equilibration time scale, more specifically, it de- 
notes the proper time point after the maximum overlap of the colliding nuclei in which 
the system reaches local thermal equilibrium. In the Bjorken model, local thermal equi- 
librium sets in at a constant proper time tq that describes a hyperbola in the Minkowski 



diagram as shown in Fig. y^. The dependence of the thermal photon spectra on tq has 
the following simple form 

oc tI (4.10) 



(Pp_i dy 

which is now derived analytically. Whenever the proper time r appeared in a model 
quantity, such as the temperature T, it was scaled by a factor of 1/tq. Thus, it is sensible 
to perform the substitution 

T = Tq W (4-11) 

in the computation of the spectra according to Eq. ( p. 281 ). By rewriting Eq. ( |3.2^ ) as 

dN „9 n 

d^p±dy 

with the integrand / denoting the production rate already integrated over the fluid cell 
rapidities y' , the result of this substitution that leads to 

= ttRIt^ / dwwf{T{w), ...) (4.13) 

d^p±dy Jwi(=Ti/To) 

can be seen very clearly as one can directly read off the ro-dependence claimed in Eq. ( [4.10| ). 
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Figure 4.7: Total Thermal Photon Spectra for tq = 0.1 fm (Lower Curves), 0.5 fm (Middle 
Curves), and 1 fm (Upper Curves) in the Phase Transition Scenario. As in the upper 
diagram of Fig. the solid lines are obtained if QGP bremsstrahlung processes are 
included and the dashed lines are obtained if QGP bremsstrahlung processes are neglected. 



The quadratical dependence ( [4.10 ) is also confirmed in Fig. 4.7, where the total thermal 
photon yields with (solid line) and without (dashed line) bremsstrahlung processes are 
displayed for three thermalization times, tq = 0.1 fm, 0.5 fm, and 1 fm. It is here 
instructive to note that each of these three values implies a different entropy, which can 
be realized by recalling that tq determines the initial volume 



Vo = 2TrR\To 



and, therefore, the entropy 



S = sqVo 

if all other parameters remain unchanged. 



(4.14) 



(4.15) 
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Figure 4.8: Lifetimes of the Collision Stages and their Dependence on the Initial Temper- 
ature Tq. The dashed, dot-dashed, and dotted lines illustrate the To-dependence of the 
lifetimes of the pure QGP, MP, and pure HHG collision stage, respectively. 

4.4.2 To - Initial Temperature 



The initial temperature Tq is the temperature of the fluid cells at the thermalization 
time tq. Before this proper time point, a temperature cannot be defined because the 
system is in its non-equilibrium phase still heading towards local thermal equilibrium. 
For the dependence of the thermal photon spectra on the initial temperature, it is not 
possible to give a relation as simple as (4.6) or ( [4.10| ), which means that the impact of 
the initial temperature must be studied on the extracted thermal photon spectra. 

The importance of the initial temperature is already anticipated by the T-dependence 
of the thermal production rates and by the To-dependence of the lifetimes of the collision 



stages. Because the evolution of the entropy density (|2.34|) provides the time scale for the 
complete hydrodynamic expansion, the cubic To-dependence of the initial entropy density. 



(4.16) 



is transferred to the lifetimes of the collision stages as can be seen in Eqs. ( |2.46| ), ( p.47| ), 
and (|2l^ ), and in Fig. |]8. 
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The above anticipation is confirmed in Figs. 4.£ and 4.1C, where the total photon 
spectra (upper diagrams) with and without bremsstrahlung processes and the photon 
spectra from the different processes in the QGP state of matter (lower diagrams) are 
illustrated for Tq = 200 MeV and Tq = 300 MeV. In the upper diagrams of these figures, 
it can be seen that an increase of the initial temperature from Tq = 200 MeV (Fig. 



up to Tq = 300 MeV (Fig. 4.10| ) leads to an increase in the total thermal photon yields by 



more than one order of magnitude and to an increase in the effect of the bremsstrahlung 
processes especially in the high-p_|_ region. It is the significantly longer lifetime of the 
thermalized collision phase and the T-dependence of the thermal production rates that 
supports high thermal photon yields in the case of high initial temperatures. The lower 
diagrams in Figs. [4.9| and 4.10| affirm the above observations concerning the contributions 



from the QGP state of matter. The QGP contribution for Tq = 300 MeV exceeds the 
one for Tq = 200 MeV by one order of magnitude at p_L = 1 GeV and by more than 
three orders of magnitude at p± = 5 GeV. A close look at the shapes of the QGP photon 
spectra is also interesting. While the Tq = 200 MeV spectra basically exhibit one slope 
that indicates a small variation in the temperature, the Tq = 300 MeV spectra exhibit 
two different slopes, where the steeper one results from the low temperature contribution 
of the MP and the flatter one corresponds to the high temperature contribution of the 
pure QGP phase. If one could really neglect transverse expansion, the occurrence of two 
different slopes in the thermal photon spectra would be a probable signature for a long 
lived MP, and thus, for the deconfinement phase transition. However, by taking into 
account also transverse expansion, the Doppler effect can mimic high temperatures and 
cause a flat slope p3, p6|. 



It is also instructive to look at Figs. 4.11 and 4.12 that present for Tq = 200 MeV and 



Tq = 300 MeV, respectively, the contributions from the two states of matter with and 
without bremsstrahlung processes taken into account. In the Tq = 200 MeV illustration 



(Fig. 4.11 ), one sees a HHG contribution that completely dominates the thermal photon 
spectrum if bremsstrahlung processes in the QGP are neglected (lower diagram). This 
changes significantly by including the production rates from the two-loop calculations 



(upper diagram). A similar behavior is observed for Tq = 300 MeV (Fig. 4.12| ), where 
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the inclusion of the QGP bremsstrahlung processes extends the region, in which the QGP 
thermal photon contribution dominates, from p± > 4: GeV to > 2 GeV. 

The comparison of the photon spectra of the different collision phases for Tq = 
200 MeV and Tq = 300 MeV manifests our above statement concerning the influence 
of the collision phase lifetime and the temperature range on the thermal photon spec- 
tra: a small lifetime of the pure QGP phase at relatively low temperatures results in 
small thermal photon yields from this phase (Fig. 4.13D , while longer lifetimes at higher 



temperatures particularly benefit the population of the high-p^ region (Fig. 4.14| ) 
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Figure 4.9: Total Thermal Photon Spectra (Upper Diagram) and Thermal Photon Spectra 
from the QGP State of Matter (Lower Diagram) in the Phase Transition Scenario with 
To = 200 MeV. Same as Fig. |3 but for Tq = 200 MeV. 
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Figure 4.10: Total Thermal Photon Spectra (Upper Diagram) and Thermal Photon Spec- 
tra from the QGP State of Matter (Lower Diagram) in the Phase Transition Scenario 
with To = 300 MeV. Same as Fig. O but for Tq = 300 MeV. 
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Figure 4.11: Thermal Photon Spectra from the QGP and the HHG State of Matter in the 



Phase Transition Scenario with Tq = 200 MeV. Same as Fig. 4.2 but for Tq = 200 MeV. 
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Figure 4.12: Thermal Photon Spectra from the QGP and the HHG State of Matter in the 



Phase Transition Scenario with Tq = 300 MeV. Same as Fig. |4j but for Tq = 300 MeV. 
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Figure 4.13: Thermal Photon Spectra from the Different CoUision Stages in the Phase 



Transition Scenario with Tq = 200 MeV. Same as Fig. 4.2 but for Tq = 200 MeV. 
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Figure 4.14: Thermal Photon Spectra from the Different Colhsion Stages in the Phase 
Transition Scenario with Tq = 300 MeV. Same as Fig. U but for Tq = 300 MeV. 
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4.5 Tc - Transition Temperature 



The transition temperature characterizes the critical point in which quarks and gluons 
become confined, or, vice versa, in which hadrons become deconfined. It should depend 
only on the nature of the strong interaction. Experimentally, the search for the transition 
temperature in ultra-relativistic heavy ion experiments is the search for discontinuities 
indicating the existence of the QGP. According to the simple model discussed in Chap. ^ 
and Eq. ( |2.33| ), particularly, a different approach in the spirit of the bag model is also 
available, where values of the bag constant B accessible through hadron spectroscopy 
can be used to infer Tc- The most sophisticated theoretical tool for the investigation of 
the transition temperature is lattice QCD, which has been applied to pure gluon theory 
[Nf = 0) and also to theories including dynamical quarks {Nf = 4, 2, 2 + 1). Depending 
upon the number of flavors embedded in the lattice QCD calculations, different values 
for the transition temperature have been found |l^ and we restrict ourselves to the 
transition temperature range 

150 MeV ~ Tc ~ 200 MeV (4.17) 
that contains most of the predicted values. 

For the dependence of the thermal photon spectra on the transition temperature Tc, 
again no simple formula can be given and the thermal photon spectra extracted within 
the simple model must be examined. Before proceeding in this way, it is instructive to 
investigate the Tc-dependence of the collision phase lifetimes Ar^, Atc, and Ar/j. There 
is no influence of Tc on the total lifetime 

At = At, + Ar^ + Ar,, (4.18) 



of the local thermal equilibrium phase since the freeze-out time given in Eq. (Ill) IS 



independent of T^- Instead, the transition temperature toggles only the distribution of 
the total lifetime on the different collision stages, pure QGP phase, MP, and pure HHG 



phase, as is illustrated in Fig. |4.15| . For low values of the transition temperature, e.g., 
Tc = 160 MeV, the lifetimes of the pure QGP phase and the MP are relatively large and 
the lifetime of the pure HHG phase is relatively short. This changes by going to high 
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Figure 4.15: Lifetimes of the Collision Stages and their Dependence on the Transition 
Temperature Tc. The dashed, dot-dashed, and dotted lines illustrate the Tc-dependence 
of the lifetimes of the pure QGP, MP, and pure HHG collision stage, respectively. By 
summing all three lifetimes in thought, it can be seen that the total collision lifetime is 
independent of Tc. 



transition temperatures, e.g., Tc = 200 MeV, where the lifetimes of the pure QGP phase 
and the MP are relatively short and the lifetime of the pure HHG phase is relatively long. 
This behavior of the lifetimes predicts the behavior of the thermal photon spectra under 
the variation of Tc'. contributions from the QGP state of matter will be higher for low 
values of Tc, while contributions from the HHG state of matter will be higher for high 



values of Tc. The thermal photon spectra shown in Figs. 4.16 and 4.17 affirm the above 
expectation. In Fig. [4.16| , one can compare the spectra from the two states of matter for 
Tc = 160 MeV (upper diagram) with those for Tc = 200 MeV (lower diagram) and see that 
the contribution from the QGP state of matter is dominant for Tc = 160 MeV, while that 
from the HHG state of matter is dominant for Tc = 200 MeV. This results in the QGP 
bremsstrahlung effect being more pronounced for Tc = 160 MeV than for Tc = 200 MeV 



as is shown in Fig. 4.17 . The total thermal photon yield displayed in both figures is also 
affected by going from Tc = 160 MeV to Tc = 200 MeV. It increases by about a factor of 
five due to a higher mean temperature in the scenario with Tc = 200 MeV. 
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Figure 4.16: Thermal Photon Spectra from the QGP and the HHG State of Matter in 
the Phase Transition Scenario for Tc = 160 MeV (Upper Diagram) and Tc = 200 MeV 



(Lower Diagram). Same as Fig. gj but for Tc = 160 MeV and = 200 MeV with QGP 
bremsstrahlung processes included in both diagrams. 
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Figure 4.17: Total Thermal Photon Spectra in the Phase Transition Scenario for Tc = 
160 MeV (Upper Diagram) and Tc = 200 MeV (Lower Diagram). Same as the upper 
diagram in Fig. |1| but for Tc = 160 MeV and Tc = 200 MeV. 
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4.6 Tf - Freeze-Out Temperature 



The freeze-out temperature Tf defines the proper time point at which the system be- 
comes dilute. At later times, the strongly interacting matter cannot be considered fluid- 
like, which means a hydrodynamical description is not justified. Instead, the picture of 
hadrons heading as free particles towards the detector is a more appropriate one for the 
phase beyond freeze-out. The definition of the freeze-out temperature is a matter of taste 



since there is no rigid criterion |24]. A common approach is however Landau's mean free 

mjp 



path consideration: when the mean free path of the hadrons A^j- r is of the order of the 



transverse size of the fireball, 



>^tfp = — ^Ra (4.19) 



with the number density Uh and the total cross section a, the system cannot be treated as 
fluid-like anymore. By applying this criterion, a freeze-out (hyper-)surface for each hadron 
species present in the strongly interacting continuum can be defined. An alternative, more 
involved way to determine Tf is to combine experimental hadron spectra with theoretical 
flow investigations. With the above means, Tf values can be found in the range ||5^ 

100 MeV ^ Tf ^ 160 MeV, (4.20) 

which is also the region over which we vary the freeze-out temperature Tf. 

As in the preceding sections, the look on the lifetimes of the collision stages gives a first 
insight into the influence of the investigated parameter. For the freeze-out temperature Tf, 
only the lifetime of the pure HHG phase expressed in Eq. ( ^.48 ) shows a dependence in 



a way that it increases with decreasing Tf. A lower value of Tf consequently leads to a 
higher thermal photon yield from the pure HHG collision phase. However, because the 
temperature at the end of the pure HHG collision phase is relatively low, the enhancement 



from the longer lived pure HHG collision phase should be relatively small. Figure 4.18 
confirms this expectation on the thermal photon yields obtained for Tf = 100 MeV 
(upper diagram) and Tf = 160 MeV (lower diagram). The marginal difference in the 
yields is only slightly visible in the low-p_|_ region. Thus, for the thermal photon spectra 
in the considered transverse momentum range, 1 GeV < p_L < 5 GeV, the infiuence of 
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Figure 4.18: Total Thermal Photon Spectra in the Phase Transition Scenario for Tf = 
100 MeV (Upper Diagram) and Tj = 160 MeV (Lower Diagram). Same as the upper 
diagram in Fig. gl| but for Tf = 100 MeV and Tf = 160 MeV. 
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the parameter Tf is weak and negligible. This fact will change if transverse expansion 



is included [24, 26 1, which should be done systematically in a future extension of this 
investigation. 



80 



Chapter 5 



Experimental Data and 
Comparison with other Works 



As mentioned in the introduction, the examination of experimental data is the most im- 
portant element of a systematic theoretical investigation. Unfortunately, the only reliable 
data on direct photon production in ultra-relativistic heavy ion reactions are the upper 
limits from fixed target 200 A ■ GeV S + Au collisions at the CERN SPS presented by the 



WA80 collaboration |16, 17|. From the successor experiment WA98, only very preliminary 
data has been presented which showed direct photon production in fixed target 158 ^-GeV 
Ph + Pb collisions at the CERN SPS. The final WA98 results are eagerly awaited and ex- 
pected to be published in the near future. In the meanwhile, the WA98 detector has been 
disassembled and the electromagnetic calorimeter LEDA (LEadglass Detector Array) has 
been shipped to the BNL RHIC, where it will be used for the measurement of the direct 
photon spectrum in the PHENIX experiment ||l^. The higher temperatures and the large 
fireball volumes that will be achieved at RHIC (due to the higher center-of-mass energy 
of y/s = 200 j4-GeV and the employed '^'^'^Au projectiles) will support higher multiplicities 
and higher thermal photon yields. RHIC is presently in its testing phase with first data 
to be taken in fall 1999. This means final direct photon data from RHIC will be available 
in about five years. Even higher energies = 5500 A -GeV) will be reached in collisions 
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of "^^^Ph nuclei at the CERN LHC that is in its construction phase and will be completed 
in the year 2005. For this machine, the dedicated heavy ion experiment ALICE (A Large 



Ion Collider Experiment) [19|] is planned. As illustrated in Fig. it will be equipped 
with the finely segmented photon spectrometer PHOS which will record the ALICE direct 
photon data. This data might be the decisive one, however, it will take about a decade 
from now for it to be analyzed and published. 

Under the above circumstances, it is clear that presently only the WA80 data can 
be inspected, while patience is required for the WA98, PHENIX, and ALICE data. We 
compare our results not only with the WA80 upper limits, but also with other theoretical 
works. There are many studies within one-fluid and even three-fluid hydro dynamical 
models that simulate also transverse expansion and employ more realistic EOS's for the 
HHG. In parallel to this thesis, the QGP bremsstrahlung processes have been considered 
in such a one-fluid 2+1 hydro dynamical model ||5^, |6^ and we will exclusively confront 
our results with the ones obtained in |^, ^] . Other investigations that do not account for 
QGP bremsstrahlung processes are also briefly reviewed. They provide hints on interesting 
aspects that should be examined in a future extension of this investigation. 



5.1 WA80 Upper Limits 



The measurement of the direct photon yields in 200 A ■ GeV S + Au collisions at the 
CERN SPS was one of the main goals in the experiment WA80. Facing a high multiplicity 
environment {dN/ dy ~ 200) and the difficulty of large backgrounds from the Dalitz decays 
TT^ and rj — > 77, the WA80 collaboration performed the analysis not event-by- 

event but on a statistical basis. In principle, the direct photon yield 7'^*'' is obtained by 
identifying the background photon yield '^^^sd subtracting it from the total observed 
photon yield 7"^^, 

^dir ^ ^ohs _ ^bkgd_ ^5 -^^ 
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Figure 5.1: The ALICE Detector. For the measurement of the direct photon spectrum, 
the photon spectrometer PHOS will be used. It is planned as a state-of-the-art electro- 
magnetic calorimeter with very high granularity. The extreme segmentation is necessary 
to handle the extreme multiplicities {dN/dy > 2500) expected due to the prospected 
position in the mid-rapidity region. 
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However, the WA80 analysis followed an alternative approach less sensitive to systematic 
error in which the photon yields were normalized with the vr*^ yield 

^ \dir / ^ \obs / 'y \ ''^9'^ 



This is reasonable since the majority of photons detected at a given p± originates from 
the decay of vr^'s at nearly the same p±. A different measure for the direct photon yield 
is also the fraction {-y / tt^)°^^ / {-y / TT^)^^9d ^j^at indicates direct photons if larger than one. 

The above expression exhibits the identification of the background photons as the 
crucial element of the direct photon analysis. Because 98% of the background photons are 
originating from the Dalitz decays vr'' — > 77 and r] 77, it is clear that the identification 
of background photons coincides with the reconstruction of vr'^ and r] mesons. Equipped 
with a lead glass photon spectrometer of fine granularity and good energy resolution, the 



WA80 detector shown in Fig. 5.2 was well suited for this task. Using the invariant mass 
method 

M^^ = {E^, +E^,)^ + {p^, + P72 )^ (5-3) 

vr*^ and r] mesons were reconstructed up to transverse momenta of several GeV. The main 
uncertainties in this method, which in fact control the systematic errors, arise from shower 
overlap due to the high photon multiplicity and from combinatorial background. The 
combinatorial background, however, could be determined with the event mixing method 
in which the invariant mass method is applied on photon pairs from artificial uncorrelated 
events. 

The final result from the WA80 analysis of central collisions is an average direct photon 
excess over background sources of 5.0% ± 0.8%(stat.) ± 5.8%(syst.) over the transverse 
momentum range 0.5 GeV < p± < 2.5 GeV. With the measured direct photon yields and 
the statistical and systematical errors obtained for central collisions, also upper limits at 
the 90% confidence level have been calculated for the direct photon yield in each p±-hm 



up to about 3.5 GeV [16|. These upper limits are the most interesting experimental results 
with regard to this investigation and we examine them within our approach to the thermal 
photon yields. 
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Figure 5.2: The WA80 Detector. Most important in the analysis of the direct photon 
production has been the photon spectrometer. It was a lead glass calorimeter of fine 
granularity and good energy resolution that covered the laboratory frame pseudorapidity 



range 2.1 < rj < 2.9 (corresponding to the mid-rapidity range in the center-of-mass frame) 
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In Chap. ^, the initial temperature Tq has been found as a model-parameter with a 
very strong effect on the thermal photon spectra. We therefore use the WA80 upper limits 
to extract an upper limit for the initial temperature Tq in the phase transition scenario 
with the remaining parameters set as follows 

Vnud = 3.0 {V^ = 20A-GeY), 

A = 32 i^^S-^^^Au), 

Qq = 37 (two-flavored QGP), 

Qh = (ideal massless pion gas), 

To = 1 fm, 

Tc = 170 MeV, 

Tf = 150 MeV. (5.4) 

Because the focus of this thesis is on the effects of the QGP bremsstrahlung processes, two 
maximum values for the initial temperature Tg"°^ are obtained such that the computed 



thermal photon yield does not exceed the WA80 upper limits. In Fig. 5.3, one sees that 
with QGP bremsstrahlung processes (solid line) a maximum initial temperature of T™"^ = 
185 MeV (upper diagram) is found, while without QGP bremsstrahlung processes (dashed 
line) a maximum initial temperature of T^""^ = 200 MeV (lower diagram) is found. By 
lowering the transition temperature Tc, these maxima can be extended to higher values 
since a lower transition temperature reduces the total thermal photon yield. (See Sec. |4.5| .) 
We illustrate the effect of a 10 MeV lower transition temperature of 

Tc = 160 MeV, (5.5) 



in Fig. 5.4, where the maximum initial temperatures T^""^ = 195 MeV (upper diagram) 
and 7g""^' = 215 MeV (lower diagram) are identified with (solid line) and without (dashed 
line) QGP bremsstrahlung contributions, respectively. The difference between the max- 
ima with and without QGP bremsstrahlung processes has grown because the QGP phase 
becomes more important with lower values of Tc and higher values of Tq. In summary, the 
inclusion of the QGP bremsstrahlung processes reduces the possible range of initial tem- 
peratures with regard to the WA80 upper limits down to lower values by about 15 MeV 
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for high Tc's and about 20 MeV for low Tc's. There is basically no effect of the QGP 
bremsstrahlung processes on the shape of the thermal photon spectrum which is a conse- 
quence of a relatively small QGP contribution in the above scenarios. For higher initial 

temperatures as expected at RHIC and LHC, the difference in the maximum initial tem- 
peratures and also in the shape of the spectrum due to the bremsstrahlung effects will be 
more pronounced. 

Additional remarks should be made on the comparison of the WA80 upper limits 
with the computed thermal photon spectra. Since the WA80 upper limits have been 
obtained from the measured direct photon yield and the corresponding statistical and 
systematical errors, not too much attention should be paid to the shape indicated by the 
upper limits. For example, the values of the upper limits in the region p±>2 GeV, which 
are high in comparison with the extracted thermal photon spectra, can be attributed to 
large statistical (less multiplicity) and large systematical (shower overlap) errors. One 
can, of course, also suspect a flattening of the direct photon spectrum due to prompt 
photons and transverse expansion, but we think that such statements inferred from the 
WA80 upper limits cannot be made without having a favor for speculation. On the other 
hand, it is interesting that the WA80 upper limits provide space for such phenomena. 
Another deviation of the shape of the thermal photon spectrum from the shape suggested 
by the upper limits is observed for p_L < 1 GeV. Although a lowering of the freeze-out 
temperature Tf would slightly reduce the gap, one should remember that the implemented 
rates have been calculated in the asymptotic regime S> T, which means that only the 
thermal photon spectra for p±>l GeV should be seriously considered. 
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Figure 5.3: WA80 Upper Limits and the Total Thermal Photon Spectra in the Phase 
Transition Scenarios {Tc = 170 MeV) with To = 185 MeV (Upper Diagram) and To = 
200 MeV (Lower Diagram). The WA80 upper limits at the 90% confidence level for the 
direct photon yield are represented in the data points, while the solid and dashed lines 
indicate the computed total thermal photon yields with and without QGP bremsstrahlung 
contributions, respectively. 
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Figure 5.4: WA80 Upper Limits and the Total Thermal Photon Spectra in the Phase 
Transition Scenarios (Tc = 160 MeV) with Tq = 195 MeV (Upper Diagram) and Tq = 
215 MeV (Lower Diagram). As in Fig. |5.3| , the data points represent the WA80 upper 
limits and the solid and dashed lines represent the computed total thermal photon yields 
with and without QGP bremsstrahlung contributions, respectively. 
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5.2 Comparison with other Works 



In parallel to this thesis, an older investigation of thermal photon production in 200 A-GeV 
S + Au collisions has been modified IgH] by implementing the QGP bremsstrahlung 



processes and a HHG EOS that included every hadron listed in the particle data table |61] 



up to a mass of 2.5 GeV. The employed model describes besides a longitudinal also a 



cylindrically symmetric transverse expansion, where the SHASTA algorithm [62| was used 
to solve the 2+1 hydrodynamic equations with the EOS. This EOS has been constructed 
as is outlined in Chap. ^, but with the mentioned much richer HHG EOS in the modified 
version. The transition and the freeze-out temperatures were set to Tc = 160 MeV 
and Tf = 120 MeV and with Eq. ([4. j)P] and dN/dy = 225, an initial temperature of 
To = 203 MeV was obtained by assuming tq = 1 fm. For the thermal photon production 
from the QGP state of matter, the rates were implemented as in this thesis, while a 
different approach was followed for the photon emission from the HHG state of matter: 
the results of Kapusta et al. [21| and Xiong et al. [27| were applied directly instead of 



using a parameterization. 

We now extract thermal photon spectra in the phase transition scenario with the 
following parameters 

Vnud = 3.0 (Vi = 20^-GeV), 
^ = 32 i^^S-^^^Au), 
Qq = 37 (two-flavored QGP), 
Qh = (ideal massless pion gas), 

To = 1 fm, 
Tq = 203 MeV, 
Tc = 160 MeV, 

Tf = 120 MeV. (5.6) 
and compare the obtained results shown in Fig. |5.5| with the results of [30| shown in 



'^For the projectile radius, the slightly different phenomenological expression Ra = 1.2 fm A^^^ was 
used. 
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Fig. 5.6. Concentrating on the region in which the implemented rates are vahd, p± > 



1 GeV, one finds basicahy identical thermal photon yields from the QGP state of matter 



(dashed lines in Fig. and dot-dashed lines in Fig. |5.6| ) and slightly different thermal 
photon yields from the HHG state of matter (dotted lines in Fig. |5.5| and dashed lines 
in Fig. ^^). The slight differences in the yields from the HHG state of matter can be 
attributed to the different implemented production rates, the different HHG EOS, and/or 
the transverse behavior of the fireball. However, one sees clearly no strong effects due 
to the transverse expansion and the richer HHG EOS, which confirms the validity of the 
simple, well understood model in the considered phase transition scenario. 



In [£0[, also the no phase transition scenario is examined and in comparison to their 
earlier work |37], in which a different HHG EOS was employed, much lower thermal 



photon yields are obtained. As has already be seen in |^6[], this is a consequence of the 
much richer HHG EOS that alters the yields only slightly in the phase transition scenario 
but significantly in the no phase transition scenario. The underlying mechanism is here, as 



discussed in Sec. |4.1.2| , the derivation of the initial temperature from the initial entropy 
density, in which also the EOS of the initial state of matter enters. For example, for 
the high number of effective degrees of freedom present in the rich HHG EOS, a lower 
To-value is obtained than for the small number of effective degrees of freedom present 
in the massless ideal pion gas. Thus, the choice of the HHG EOS determines the initial 
temperature in the no phase transition scenario, which has a significant influence on the 
thermal photon spectra. 



The approach used in |]60[ was also applied in ||59|] , where thermal photon spectra are 
examined for central Pb + Pb collisions at SPS, RHIC, and LHC under the assumption 
of a i/iree-flavored QGP state of matter (Nf = 3) and a lower freeze-out temperature of 
Tf = 100 MeV. Table shows the initial temperatures Tq that were obtained^ for the 



displayed assumptions of the thermalization time tq and the estimated dN/dy entries [SC]. 
To compute comparable thermal photon spectra in our approach, modifications in the 
QGP thermal photon production rates are necessary due to the assumed three- flavored 



^The initial temperatures were again calculated from Eq. (4.8) but this time with at = aq = 47.5 7r^/90 
according to Nf = 3. 
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Figure 5.5: WA80 Upper Limits and the Thermal Photon Spectra from the QGP and the 
HHG State of Matter. The thermal photon yields from the QGP and the HHG state of 
matter are represented in the dashed and the dotted lines, respectively, while the sum of 
both is the total thermal photon yield represented in the solid line. QGP bremsstrahlung 
processes are included in the upper plot and neglected in the lower plot. The WA80 upper 
limits are indicated in the data points. 
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Figure 5.6: WA80 Upper Limits and the Thermal Photon Spectra from the QGP and the 

HHG State of Matter Extracted by Srivastava and Sinha [^] . The thermal photon yields 

from the QGP and the HHG state of matter are represented in the dot-dashed (QM) and 

the dashed (HM) lines, respectively, while the sum of both is the total thermal photon 

yield represented in the solid (Sum) line. QGP bremsstrahlung processes are included in 

the upper plot and neglected in the lower plot. The WA80 upper limits are indicated in 

the data points with the arrows pointing downwards. A parton cascade model result for 
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direct photon production is illustrated in the histogram (PCM). 
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Table 5.1: Accelerator Specific Quantities Used by Srivastava for the Extraction of Ther- 



mal Photon Spectra |5^]. The dN/dy entries were estimated as proposed in [56| and 
the To assumptions were made in light of parton cascade model results |5^, Equa- 
tion (O) with = = 47.5 7r^/90 was used to compute the initial temperatures from 



the corresponding values of dN/dy and tq. 

QGP: the sum of the squared electrical charge numbers includes now also the strange 
quark and the difference of the A'^^-dependent integrals Jt and Jl becomes |5^ 

{Jt - Jl) = 9.32. (5.7) 



With these modifications, the accelerator specific values listed in Tab. 5.1, and the re- 
maining model parameters set to 



9q 

Qh 



47.5 (three-fiavored QGP), 
3 (ideal massless pion gas). 
160 MeV, 
100 MeV. 



(5.8) 



we receive the thermal photon spectra shown in the lower halfs of Figs. 5^, 5^, and 5.E . 
The upper halfs of these figures present the results of |59] to allow a direct comparison. 
At SPS energies, the thermal spectra from the QGP state of matter are again basically 
identical for p± > 1 GeV, while the thermal spectra from the HHG state of matter 
show some differences especially in the high p_|_-range. This difference becomes even more 
pronounced at RHIC and LHC energies and is a consequence of the different transverse 
behavior of the fireball: the transverse expansion modeled in |5^ causes a Doppler effect 
that mimics a higher temperature in the HHG. Also for the QGP state of matter, a slight 
blueshift in the spectra can be seen at RHIC and LHC energies. 
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Figure 5.7: Thermal Photon Spectra from Central Pb + Pb Collisions at SPS Energies. 
The upper diagram is taken from (Upper Diagram), while the lower diagram has 
been calculated in our approach. In both diagrams, the contributions from the HHG 
("Had Mat" in the upper diagram and ^Hotalh" (dotted line) in the lower diagram), 
the QGP with bremsstrahlung processes {^^QMj^urenche et at" in the upper diagram and 
^^totalq" (solid line) in the lower diagram), and the QGP without bremsstrahlung processes 
( ^^QMxapusta et al^ in the Upper diagram and by "1 — /oop" (dashed line) in the lower 
diagram) are illustrated. In addition, the 9gper diagram shows also an expectation for 
the prompt photon yield which is labeled by "QCD" . 



Pb + Pb @ RHIC; Hydrodynamics 




P± [GeV] 



Photon Spectrum (Pb+Pb @ RHIC) 



I 

> 

o 

Si 

H 

4" 



10-^ - 



10- 



10- 



1 1 


1 


1 

total, 






1 — loop 






total/, 


— >e 

N. ^ 


To 


= 0.5 fm 




To 


= 310 MeV - 






= 160 MeV - 






= 100 MeV- 








- y = o 






1 1 


1 





2 3 

P± [GeV] 



Figure 5.8: Thermal Photon Spectra from Central Pb + Pb Collisions at RHIC Energies. 



Same as Fig. but for RHIC energies. 
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Figure 5.9: Thermal Photon Spectra from Central Pb + Pb Collisions at LHC Energies. 
Same as Fig. but for LHC energies. 
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There are many other studies of thermal photon production in ultra-relativistic heavy 
ion colhsions besides the two cited above. However, they were performed without taking 



the QGP bremsstrahlung processes into account. In |24] and |26|, the 2+1 hydrodynam- 
ical model based on the SHASTA algorithm was used to study the transverse expansion 
and the effects of the different HHG EOS's on the thermal photon spectra, where the 



initial conditions were obtained as in |£S|, 6C]. A different approach to the initial con 



ditions was followed in Sollfrank et al. extracted also hadron spectra and fitted 

them with the choice of the initial conditions to the experimentally measured spectra. 
Then, these initial conditions were used in the calculation of the thermal photon spectra. 
In addition, the influence of different HHG EOS's including finite baryon density was 
checked. The underlying model of this study was again the 2+1 hydro dynamical calcula- 
tion performed with the SHASTA algorithm. An alternative algorithm which solves 3+1 
hydro dynamical equations with the EOS is HYLANDER that was applied in |23]. 



Within this simulation, no initial conditions need to be specified because the initial com- 
pression (or collision) stage is also modeled. Also in three-fluid hydrodynamical models, 
the compression stage is simulated and no initial conditions must be specified. Three-fluid 
hydrodynamics describes both of the colliding nuclei and the central region of secondary 
particles with separate fluids. In this spirit, thermal photons have been used as a measure 
for the rapidity dependence of the temperature []6^ ]. 

Most of the above investigations have been performed for the experiments at the 
CERN SPS, WA80 and WA98, in which the assumption of a chemically equilibrated 
plasma should be a good approximation, while the assumption of a vanishing baryon 
density should be only a crude approximation. A complementary picture seems to arise 
for RHIC and LHC, where one should face a central region of zero baryon density but 
no chemically equilibrated plasma. Therefore, methods have been developed to describe 
a chemically non-equilibrated QGP |66], which have already been applied in extractions 



of thermal photon spectra at RHIC and LHC [48, |5^. An alternative equilibration 
scenario is the one of the increasingly strongly interacting parton plasma (ISIPP) in which 
photon production has also been studied [p^] . 
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Chapter 6 



Conclusion 



A systematic investigation of hard thermal photon spectra from ultra-relativistic heavy 
ion reactions was presented with emphasis on the effects of brcmsstrahlung processes in 
the quark-gluon plasma (QGP). The thermal photon spectra were computed from ther- 
mal photon production rates within a simple, well understood model for central ultra- 
relativistic heavy ion collisions in the transparent energy region. This model, which is 
based on the Bjorken model, describes the fireball evolution in one-fluid 1-1-1 hydrody- 
namics as an only longitudinally expanding tube in which entropy is conserved (adiabatic 
expansion). With the implemented equation of state (EOS) that is an ingredient in 
every hydrodynamic calculation, a phase transition and a no phase transition scenario 
were simulated for the central region with vanishing baryon density. In the phase tran- 
sition scenario, a first-order phase transition from QGP to hot-hadronic gas (HHG) was 
implemented, where the QGP and the HHG were modeled respectively by an ideal mass- 
less parton gas and an ideal massless pion gas. The matching at the critical boundary 
was achieved with the Gibbs criteria and Maxwell construction. The parameters that 
determined the evolution in this scenario were the thermalization time tq, the initial tem- 
perature To, the transition temperature Tc, and the freeze-out temperature Tf. In the no 
phase transition scenario, only the HHG state of matter was assumed and modeled using 
the ideal massless pion gas. Here, the parameters tq, Tq, and Tf described the fireball 
evolution completely. Two additional parameters had to be specified in the extraction of 
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the thermal photon spectra, the mass number A of the projectile determining the trans- 
verse size of the fireball and the projectile rapidity ynud determining the longitudinal 
limits of the fireball. Finally, the thermal photon production rates entered as the most 
critical input in the computation of the thermal photon spectra. For the QGP, the 1-loop 
production rates, Compton scattering and gg-annihilation, and the recently derived 2-loop 
production rates describing "ordinary" bremsstrahlung (bremss) and gg-annihilation with 
an additional scattering in the medium (qq—aws) were considered, while for the HHG, the 
rate describing the processes irp — ^ ttj and 7rp — ^ ai — ^ ir^ was implemented. Our prime 
interest here and throughout the complete thesis was on the 2-loop production rates in 
the QGP that contain the bremsstrahlung processes. In fact, the astonishing result that 
these rates are in the same order of the coupling constants as the 1-loop rates triggered the 
research for this thesis. The main issue was then on the impact of these bremsstrahlung 
processes on the thermal photon spectrum. It was attacked in a systematic investigation, 
in which the influence of every model parameter was checked. Already the photon spec- 
tra at fixed temperatures amplified the interest in the bremsstrahlung processes because 
the qq — aws process appeared as the dominant thermal photon source over the 1-loop 
processes. 

A first look at the thermal photon spectra additionally affirmed this observation. 
In the phase transition sample scenario with the parameters set to tq = 1 fm, Tq = 
250 MeV, Tc = 170 MeV, Tf = 150 MeV, ynud = 8.6, and A = 208, a significant 
enhancement of the thermal photon yield was found in the high-p^ region, p± > 3 GeV, 
which could be traced back to an enhancement of about one order in magnitude over 
the complete considered p^-range in the yield from the QGP phase due to the two- 
loop processes taken into account. This enhancement of the yields from QGP state of 
matter was found in the subsequent systematic investigation to be independent of specific 
parameter settings. Thus, the strength of this two-loop effect coincides with the strength 
of the QGP contribution, which goes up for high initial temperatures Tq and low transition 
temperatures Tc. In summary, the following dependences were found by considering each 
model parameter, mass number A of projectile, projectile rapidity ynud^ thermalization 
time To, initial temperature Tq, transition temperature Tc, and freeze-out temperature Tf, 
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separately, with the above phase transition sample scenario serving as the basis for the 
investigation. The spectra turned out to be proportional to J^l"^ and Tq and nearly 
independent of ynud and Tf. For Tq and the dependence could not be expressed 
in terms of a simple proportionality. Instead, the different collision phase lifetimes and 
the computed spectra were closely examined. The strongest dependence of the thermal 
photon yields with respect to their magnitude and their slope was found for the initial 
temperature Tq. Because the collision phase lifetimes show basically a cubic dependence 
on To and the photon spectra show basically a quadratic dependence on T, this is clearly 
understood. In fact, Tq is the crucial parameter for the total thermal photon yields. The 
dependence on is not so significant but also interesting: for low values of Tc, the pure 
QGP phase is more dominant than for high values of T^. This was directly observed in 
the lifetimes of the collision phases, where a low value of supports a relatively long 
lived QGP phase and relatively short lived mixed and pure HHG phases. An influence on 
the total spectrum has also been observed: higher values of result in higher thermal 
photon yields. This behavior could be traced back to a mean temperature that is higher 
for high values of Tc than for low values of Tc- 

The comparison of the thermal photon spectra from the phase transition scenario with 
those from the no phase transition scenario, which should exhibit the quality of thermal 
photons as a potential signature for the QGP formation, was only briefly considered 
since the ideal massless pion gas was found inappropriate to describe the purely hadronic 
scenario. Thus, this important comparison is recommended for a future extension of this 
work, in which a more realistic EOS should be used to describe the HHG state of matter. 

The look at experimental data was unfortunately limited by the availability of di- 
rect photon data. So far, only upper limits for the direct photon production have been 
extracted from fixed target 200 A ■ GeV S + Au data at the CERN SPS. They were 
used to identify maximum values for the initial temperature, where the value found with 

bremsstrahlung processes was 15 MeV (for high Tc's) to 20 McV (for low T^i's) below the 
value found without bremsstrahlung processes. These are no severe restrictions on the 
model parameters, but they illustrate the effect of the bremsstrahlung processes in terms 
of a difference in the initial temperature. For more severe restrictions on the parameters, 
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one must wait for the WA98, PHENIX, and ALICE data which should be more exphcit 
because of the advanced experimental instrumentation. 

While waiting for the WA98 results on direct photon production, other theoretical 
works were considered, where the comparison with the investigation of Srivastava and 
Sinha was particularly informative. In parallel to this thesis, Srivastava and Sinha in- 
spected the 2-loop bremsstrahlung processes in a hydrodynamical model that describes 
also transverse expansion (2+1 hydrodynamics) and implements a more realistic HHG 
EOS. Interestingly, their results for the phase transition scenario corresponding to the 
200 A ■ GeV S + Au collisions at the CERN SPS were almost identical with ours. This 
means that in the considered scenario, the simple 1+1 hydrodynamical model is as com- 
petent as the more sophisticated 2+1 hydrodynamical model. For different scenarios 



describing Pb + Pb collisions at SPS, RHIC, and LHC energies, another recent study |5£] 
was found on the preprint server that used again the 2+1 hydrodynamical model with the 
more realistic HHG EOS to explore the effects of the bremsstrahlung processes. In com- 
parison to this study, some significant differences were observed especially for the HHG 
photon contribution in the high-p_|_ region, which were traced back to the different trans- 
verse behavior. These differences indicate that the domain of the late collision stage in a 
long lived fireball is sensible to transverse behavior. Other theoretical investigation, which 
did not include the QGP bremsstrahlung processes, were also briefly reviewed. They pro- 
vide hints on further interesting physics aspects that should be examined systematically 
in a future extension of this investigation. 

In prospect of the upcoming experimental data, we strongly recommend a step-by- 
step extension of this investigation, in which the very systematical approach should be 
kept as the first principle. First, more realistic EOS's for the QGP and, as already 
mentioned, for the HHG state of matter should be implemented. For the QGP, one 
could use parameterizations of lattice QCD data or take into account effective parton 
masses |^, For the HHG, finite hadron masses and a higher number of hadrons 

should be included. This will allow the important comparison of the results from the 
phase transition scenario with those from the no phase transition scenario, which has 
been beyond the scope of this thesis due to the inappropriate HHG EOS of the ideal 
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massless pion gas. Second, the EOS and also the photon production rates should be 
extended to the case of finite baryon density. In view of the expected WA98 data for 
Pb + Pb collisions at the CERN SPS, it will be important to understand the influence 
of finite baryon density. For RHIC and LHC, finite baryon density calculations will not 
be necessary since the degree of stopping is predicted much weaker than for SPS, but on 
the other hand, chemical equilibrium will not be fulfilled. Thus, the third step should 
be the extension of the hydrodynamical model and the photon production rates to the 
case of chemical non-equilibrium, which could be done by using fugacities. When a solid 
understanding of the above phenomena and their effects on the thermal photon spectra 
has been collected, we propose to include transverse expansion as the final step, where one 
could use the SHASTA algorithm (2+1), the HYLANDER algorithm (3+1) or develop a 
different new approach. Once the systematic investigation is driven up to this point for 
thermal photons, it should be easily transferable to dileptons and hadrons. 

With the heavy ion physics program underway at RHIC and LHC, the above study 
would definitively be of great interest over the next decade and could participate in the 

potential discovery of the QGP formation at RHIC or LHC. The C++ program developed 
for this thesis could serve as a solid basis for the above studies since it was written in an 
object oriented style ready for extension. 
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